ecoblender

alex.filazzola

Abstract

Shrubs facilitate the abundance and productive of annual plants in desert ecosystems.However, these shrub microhabitats favour plant species with competitive life histories. Ephedra californica is a dominant shrub in the San Joaquin desert that has been identified as a facilitator. Here, we explore the factors that limit the Ephedra californica recruitment into the San Joaquin desert including substrate, water availability, and herbivory. We also explore the role of the invasive grass Bromus madritensis on limiting establishment of E. californica. Vegetation surveys were conducted in the field during 2013 and collected seed was then used to conduct two greenhouse trials. The first explore germination and establishment techniques for E. california and the optimal substrate. The second examined E. californica establishment in present of the invasive B. madritensis responding to different water levels and herbivory. These results can have implications for land managers in the San Joaquin Valley to maintain native shrub biodiversity in the region.


## load libraries
library(tidyverse)
library(OIsurv)
library(lsmeans)
library(ggthemes)

## load data
substrate <-read.csv("data/Ephedra.substrate.csv")
recruit <-read.csv("data/Ephedra.recruitment.csv")
landscape <-read.csv("data/ephedra.landscape.csv")

##load functions
## inverse hyperbolic sine transformation 
ihs <- function(x) {
    y <- log(x + sqrt(x ^ 2 + 1))
    return(y)
}
se <- function(x, ...) sqrt(var(na.omit(x))/length(na.omit(x)))
source("functions.r")

Ephedra characteristics at the landscape

avg.size <- landscape %>% summarize(H=mean(H), D1=mean(D1), Area=mean(Area), density=mean(Shrub.density), rdm=mean(RDM.2013))
avg.size
##          H       D1     Area     density      rdm
## 1 1.332925 3.446412 9.321902 0.004318859 7.042223

Ephedra recruitment at landscape

## compare ephedra size and density across landscape
landscape[,"area.group"] <- as.numeric(cut(landscape$Log.area,20))
hist(landscape$area.group)

mean.area <- landscape %>% group_by(area.group) %>%  summarize(rdm=mean(RDM.2013), area=mean(Log.area), density=mean(Shrub.density))

ggplot(mean.area) + geom_point(aes(x=density, y=area), size=3)+ylab("log (area of shrub)")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16)) + xlab(expression("shrub density (m"^"2"*")")) +  stat_smooth(method="lm", formula= y~x,aes(x=density, y=area))

m1 <- lm(area~density, data=mean.area)
summary(m1)
## 
## Call:
## lm(formula = area ~ density, data = mean.area)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.82510 -0.26068  0.01197  0.22216  1.16732 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)     9.3959     0.9852   9.537 1.84e-08 ***
## density     -1706.0323   223.6888  -7.627 4.81e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5569 on 18 degrees of freedom
## Multiple R-squared:  0.7637, Adjusted R-squared:  0.7506 
## F-statistic: 58.17 on 1 and 18 DF,  p-value: 4.806e-07
ggplot(mean.area) + geom_point(aes(x=density, y=rdm), size=mean.area$area+1)+ylab("residual dry matter")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16)) + xlab(expression("shrub density (m"^"2"*")")) +  stat_smooth(method="lm", formula= y~exp(x),aes(x=density, y=rdm))

m2 <- lm(log(rdm)~density, data=mean.area)
summary(m2)
## 
## Call:
## lm(formula = log(rdm) ~ density, data = mean.area)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.32696 -0.09127 -0.02044  0.05849  0.46804 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.1851     0.3294   9.670 1.49e-08 ***
## density     -279.3744    74.7892  -3.735  0.00151 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.1862 on 18 degrees of freedom
## Multiple R-squared:  0.4367, Adjusted R-squared:  0.4054 
## F-statistic: 13.95 on 1 and 18 DF,  p-value: 0.001514

Ephedra optimal substrate

substrate[,"Sand"] <- as.factor(substrate$Sand)
substrate[,"census"] <- as.factor(substrate$census)

m1 <- aov(height ~ Micro * Sand,  data=subset(substrate, census==10))
summary(m2) ## nothing significant
## 
## Call:
## lm(formula = log(rdm) ~ density, data = mean.area)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.32696 -0.09127 -0.02044  0.05849  0.46804 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    3.1851     0.3294   9.670 1.49e-08 ***
## density     -279.3744    74.7892  -3.735  0.00151 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.1862 on 18 degrees of freedom
## Multiple R-squared:  0.4367, Adjusted R-squared:  0.4054 
## F-statistic: 13.95 on 1 and 18 DF,  p-value: 0.001514
## compare shade on survival of ephedra
my.surv <- Surv(as.numeric(substrate$census), substrate$survival)
fit1 <- survfit(my.surv~Micro, data=substrate)
summary(fit1)
## Call: survfit(formula = my.surv ~ Micro, data = substrate)
## 
##                 Micro=Open 
##  time n.risk n.event survival std.err lower 95% CI upper 95% CI
##     1    960      11    0.989 0.00343        0.982        0.995
##     2    880      29    0.956 0.00681        0.943        0.969
##     3    800      41    0.907 0.00987        0.888        0.927
##     4    720      42    0.854 0.01221        0.830        0.878
##     5    640      67    0.765 0.01504        0.736        0.795
##     6    481      31    0.715 0.01647        0.684        0.748
##     7    400      62    0.604 0.01901        0.568        0.643
##     8    241      31    0.527 0.02108        0.487        0.570
##     9    160      31    0.425 0.02366        0.381        0.474
##    10     80      31    0.260 0.02730        0.212        0.320
## 
##                 Micro=Shade 
##  time n.risk n.event survival std.err lower 95% CI upper 95% CI
##     1    960      33    0.966 0.00588        0.954        0.977
##     2    880      43    0.918 0.00897        0.901        0.936
##     3    800      42    0.870 0.01117        0.849        0.892
##     4    720      40    0.822 0.01290        0.797        0.848
##     5    640      34    0.778 0.01422        0.751        0.807
##     6    480      17    0.751 0.01521        0.721        0.781
##     7    400      31    0.692 0.01725        0.659        0.727
##     8    240      16    0.646 0.01958        0.609        0.686
##     9    160      15    0.586 0.02317        0.542        0.633
##    10     80      13    0.491 0.03099        0.433        0.555
par(mar=c(4.5,4.5,.5,.5))
plot(fit1, col="white", ylim=c(0.2,1), xlim=c(0.9,10.1), ylab="Estimated Survival Function", xlab="Census", cex.lab=1.5, cex.axis=1.3)
lines(1:10,summary(fit1)[[10]][1:10], col="#E69F00", lty=2) ## upper
lines(1:10,summary(fit1)[[6]][1:10], col="#E69F00", lwd=2) ## value sun
lines(1:10,summary(fit1)[[11]][1:10], col="#E69F00", lty=2) ## lower
lines(1:10,summary(fit1)[[10]][11:20], col="#56B4E9", lty=2) ## upper
lines(1:10,summary(fit1)[[6]][11:20], col="#56B4E9", lwd=2) ## value shade
lines(1:10,summary(fit1)[[11]][11:20], col="#56B4E9", lty=2) ## lower
legend(1.8, 0.35, c("Sun","Shade"), lty=1, lwd=3, col=c("#E69F00","#56B4E9"), cex=1.5)

coxph.fit <- coxph(my.surv ~ Micro, method="breslow", data=substrate)
coxph.fit 
## Call:
## coxph(formula = my.surv ~ Micro, data = substrate, method = "breslow")
## 
##               coef exp(coef) se(coef)     z       p
## MicroShade -0.2802    0.7557   0.0786 -3.56 0.00037
## 
## Likelihood ratio test=12.8  on 1 df, p=0.000342
## n= 1920, number of events= 660
## compare sand on survival of ephedra
my.surv <- Surv(as.numeric(substrate$census), substrate$survival)
fit2 <- survfit(my.surv~Sand, data=substrate)
summary(fit2)
## Call: survfit(formula = my.surv ~ Sand, data = substrate)
## 
##                 Sand=0 
##  time n.risk n.event survival std.err lower 95% CI upper 95% CI
##     1    384       7    0.982 0.00683        0.968        0.995
##     2    352      11    0.951 0.01125        0.929        0.973
##     3    320      13    0.912 0.01506        0.883        0.942
##     4    288      12    0.874 0.01799        0.840        0.910
##     5    256      14    0.827 0.02106        0.786        0.869
##     6    192       6    0.801 0.02289        0.757        0.847
##     7    160      12    0.741 0.02695        0.690        0.795
##     8     96       6    0.694 0.03120        0.636        0.758
##     9     64       6    0.629 0.03794        0.559        0.708
##    10     32       6    0.511 0.05325        0.417        0.627
## 
##                 Sand=25 
##  time n.risk n.event survival std.err lower 95% CI upper 95% CI
##     1    384      10    0.974 0.00813        0.958        0.990
##     2    352      16    0.930 0.01331        0.904        0.956
##     3    320      20    0.872 0.01772        0.838        0.907
##     4    288      18    0.817 0.02075        0.777        0.859
##     5    256      15    0.769 0.02292        0.726        0.815
##     6    192       7    0.741 0.02441        0.695        0.791
##     7    160      14    0.676 0.02776        0.624        0.733
##     8     96       7    0.627 0.03137        0.568        0.692
##     9     64       7    0.558 0.03714        0.490        0.636
##    10     32       7    0.436 0.05007        0.348        0.546
## 
##                 Sand=50 
##  time n.risk n.event survival std.err lower 95% CI upper 95% CI
##     1    384      11    0.971 0.00851        0.955        0.988
##     2    352      16    0.927 0.01350        0.901        0.954
##     3    320      19    0.872 0.01765        0.838        0.907
##     4    288      20    0.812 0.02098        0.771        0.854
##     5    256      18    0.755 0.02343        0.710        0.802
##     6    192       9    0.719 0.02512        0.672        0.770
##     7    160      17    0.643 0.02848        0.589        0.701
##     8     96       8    0.589 0.03178        0.530        0.655
##     9     64       8    0.516 0.03697        0.448        0.593
##    10     32       8    0.387 0.04823        0.303        0.494
## 
##                 Sand=75 
##  time n.risk n.event survival std.err lower 95% CI upper 95% CI
##     1    384       9    0.977 0.00772        0.962        0.992
##     2    352      18    0.927 0.01361        0.900        0.954
##     3    320      16    0.880 0.01716        0.847        0.915
##     4    288      17    0.828 0.02025        0.790        0.869
##     5    256      28    0.738 0.02422        0.692        0.787
##     6    192      13    0.688 0.02624        0.638        0.741
##     7    160      26    0.576 0.02976        0.521        0.637
##     8     96      13    0.498 0.03266        0.438        0.566
##     9     64      13    0.397 0.03612        0.332        0.474
##    10     32      13    0.236 0.04058        0.168        0.330
## 
##                 Sand=100 
##  time n.risk n.event survival std.err lower 95% CI upper 95% CI
##     1    384       7    0.982 0.00683        0.968        0.995
##     2    352      11    0.951 0.01125        0.929        0.973
##     3    320      15    0.907 0.01554        0.877        0.937
##     4    288      15    0.859 0.01891        0.823        0.897
##     5    256      26    0.772 0.02349        0.727        0.819
##     6    193      13    0.720 0.02596        0.671        0.773
##     7    160      24    0.612 0.03000        0.556        0.674
##     8     97      13    0.530 0.03351        0.468        0.600
##     9     64      12    0.431 0.03755        0.363        0.511
##    10     32      10    0.296 0.04372        0.222        0.395
test <- substrate %>% group_by(census,Sand) %>%  summarize(count=sum(survival),avg=mean(survival))
test <- data.frame(test)

par(mar=c(4.5,4.5,.5,.5))
plot(fit2, col="white", ylim=c(0.2,1), xlim=c(0.9,10.1), ylab="Estimated Survival Function", xlab="Census", cex.lab=1.5, cex.axis=1.3)
lines(1:10,summary(fit2)[[10]][1:10], col="#E69F00", lty=2) ## upper
lines(1:10,summary(fit2)[[6]][1:10], col="#E69F00", lwd=2) ## value sand 0
lines(1:10,summary(fit2)[[11]][1:10], col="#E69F00", lty=2) ## lower
lines(1:10,summary(fit2)[[10]][11:20], col="#56B4E9", lty=2) ## upper
lines(1:10,summary(fit2)[[6]][11:20], col="#56B4E9", lwd=2) ## value sand 25
lines(1:10,summary(fit2)[[11]][11:20], col="#56B4E9", lty=2) ## lower
lines(1:10,summary(fit2)[[10]][21:30], col="#009E73", lty=2) ## upper
lines(1:10,summary(fit2)[[6]][21:30], col="#009E73", lwd=2) ## value sand 50
lines(1:10,summary(fit2)[[11]][21:30], col="#009E73", lty=2) ## lower
lines(1:10,summary(fit2)[[10]][31:40], col="#CC79A7", lty=2) ## upper
lines(1:10,summary(fit2)[[6]][31:40], col="#CC79A7", lwd=2) ## value sand 75
lines(1:10,summary(fit2)[[11]][31:40], col="#CC79A7", lty=2) ## lower
lines(1:10,summary(fit2)[[10]][41:50], col="#D55E00", lty=2) ## upper
lines(1:10,summary(fit2)[[6]][41:50], col="#D55E00", lwd=2) ## value sand 100
lines(1:10,summary(fit2)[[11]][41:50], col="#D55E00", lty=2) ## lower
legend(1.7, 0.44, c("0%","25%","50%","75%","100%"), lty=1, lwd=3, col=c("#E69F00","#56B4E9","#009E73","#CC79A7","#D55E00"), cex=1.2, title="Sand")

coxph.fit <- coxph(my.surv ~ Sand, method="breslow", data=substrate)
coxph.fit 
## Call:
## coxph(formula = my.surv ~ Sand, data = substrate, method = "breslow")
## 
##          coef exp(coef) se(coef)    z       p
## Sand25  0.263     1.301    0.138 1.91  0.0563
## Sand50  0.365     1.441    0.135 2.71  0.0068
## Sand75  0.579     1.785    0.130 4.47 7.7e-06
## Sand100 0.450     1.568    0.133 3.39  0.0007
## 
## Likelihood ratio test=23.3  on 4 df, p=0.000109
## n= 1920, number of events= 660
## raw number of germinants end trials

end <- subset(substrate, census==10)

## compare microsite

end.means <- end %>% group_by(Micro) %>% summarize(eph=mean(survival), eph.se=se(survival), Height=mean(height, na.rm=T), Height.se=se(height)) %>% gather(measure, value, 2:5) %>% separate(measure, c("Eph",".se")) ## calculate metrics and put long format
## Warning: Too few values at 4 locations: 1, 2, 5, 6
end.means[is.na(end.means$.se),3] <- "mean"
end.means <- end.means %>% spread(.se, value)
end.means <- data.frame(end.means)

end.means[,"lower"] <- end.means[,"mean"]-end.means[,"se"]
end.means[,"upper"] <- end.means[,"mean"]+end.means[,"se"]

survival <- subset(end.means, Eph=="eph")

plot1 <- ggplot(survival, aes(x=Micro, y=mean)) + geom_bar(stat="identity") +
  geom_errorbar(aes(ymin=lower, ymax=upper), width=.2)+ylab("average surviving Ephedra plant")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16)) + xlab("Microsite") 

height <- subset(end.means, Eph=="Height")

plot2 <- ggplot(height, aes(x=Micro, y=mean)) + geom_bar(stat="identity") +
  geom_errorbar(aes(ymin=lower, ymax=upper), width=.2)+ylab("average height of Ephedra plant")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16)) + xlab("Microsite") 

require(gridExtra)
## Loading required package: gridExtra
## 
## Attaching package: 'gridExtra'
## The following object is masked from 'package:dplyr':
## 
##     combine
grid.arrange(plot1, plot2, ncol=2)

m1 <- glm(survival~ Sand* Micro, data=end, family="binomial")
anova(m1, test="Chisq")
## Analysis of Deviance Table
## 
## Model: binomial, link: logit
## 
## Response: survival
## 
## Terms added sequentially (first to last)
## 
## 
##            Df Deviance Resid. Df Resid. Dev Pr(>Chi)   
## NULL                         159     188.21            
## Sand        4   4.7398       155     183.47 0.315056   
## Micro       1  10.7288       154     172.75 0.001055 **
## Sand:Micro  4   8.6588       150     164.09 0.070217 . 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
lsmeans(m1, pairwise~Sand)
## NOTE: Results may be misleading due to involvement in interactions
## $lsmeans
##  Sand     lsmean          SE df   asymp.LCL    asymp.UCL
##  0    -1.4663371   0.4529108 NA   -2.354026  -0.57864819
##  25   -1.3671838   0.4643107 NA   -2.277216  -0.45715153
##  50   -8.7830342 494.5226042 NA -978.029528 960.46345952
##  75   -0.4236489   0.3831780 NA   -1.174664   0.32736619
##  100  -0.8047190   0.3872983 NA   -1.563810  -0.04562817
## 
## Results are averaged over the levels of: Micro 
## Results are given on the logit (not the response) scale. 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast    estimate          SE df z.ratio p.value
##  0 - 25   -0.09915331   0.6486236 NA  -0.153  0.9999
##  0 - 50    7.31669717 494.5228116 NA   0.015  1.0000
##  0 - 75   -1.04268814   0.5932568 NA  -1.758  0.3987
##  0 - 100  -0.66161811   0.5959263 NA  -1.110  0.8012
##  25 - 50   7.41585049 494.5228222 NA   0.015  1.0000
##  25 - 75  -0.94353482   0.6020048 NA  -1.567  0.5185
##  25 - 100 -0.56246480   0.6046358 NA  -0.930  0.8852
##  50 - 75  -8.35938531 494.5227526 NA  -0.017  1.0000
##  50 - 100 -7.97831529 494.5227559 NA  -0.016  1.0000
##  75 - 100  0.38107003   0.5448168 NA   0.699  0.9567
## 
## Results are averaged over the levels of: Micro 
## Results are given on the log odds ratio (not the response) scale. 
## P value adjustment: tukey method for comparing a family of 5 estimates
## survival
end.means <- end %>% group_by(Sand) %>% summarize(eph=mean(survival), eph.se=se(survival))
end.means <- data.frame(end.means)

end.means[,"lower"] <- end.means[,"eph"]-end.means[,"eph.se"]
end.means[,"upper"] <- end.means[,"eph"]+end.means[,"eph.se"]

ggplot(end.means, aes(x=Sand, y=eph)) + geom_bar(stat="identity", fill="#56B4E9") +
  geom_errorbar(aes(ymin=lower, ymax=upper), width=.2)+ylab("average surviving Ephedra plant")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16)) + xlab("Percentage of sand in soil") 

## height

end.means <- end %>% group_by(Sand) %>% summarize(eph=mean(height, na.rm=T), eph.se=se(height, na.rm=T))
end.means <- data.frame(end.means)

end.means[,"lower"] <- end.means[,"eph"]-end.means[,"eph.se"]
end.means[,"upper"] <- end.means[,"eph"]+end.means[,"eph.se"]

ggplot(end.means, aes(x=Sand, y=eph)) + geom_bar(stat="identity", fill="#56B4E9") +
  geom_errorbar(aes(ymin=lower, ymax=upper), width=.2)+ylab("average surviving Ephedra plant")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16)) + xlab("Percentage of sand in soil") 

## most germinants
maximum.germ <- substrate %>% group_by(Micro, Sand, Rep) %>% summarize(eph.max=max(survival))

max.eph <- maximum.germ %>% group_by(Sand) %>% summarize(eph=mean(eph.max,na.rm=T),eph.se=se(eph.max, na.rm=T)) 
max.eph <- data.frame(max.eph)

max.eph[,"lower"] <- max.eph[,"eph"]-max.eph[,"eph.se"]
max.eph[,"upper"] <- max.eph[,"eph"]+max.eph[,"eph.se"]


ggplot(max.eph, aes(x=Sand, y=eph)) + geom_point(fill="black", size=4) +
  geom_errorbar(aes(ymin=lower, ymax=upper), width=.2)+ylab("total Ephedra germinants")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16)) + xlab("Percentage of sand in soil") + ylim(0,1) + stat_smooth(method="lm", formula= y~poly(x,2),data=max.eph, aes(x=as.numeric(Sand), y=eph), lwd=1, color="black", lty=2)

m1 <- lm(eph~poly(as.numeric(Sand),2), data=max.eph)
summary(m1)
## 
## Call:
## lm(formula = eph ~ poly(as.numeric(Sand), 2), data = max.eph)
## 
## Residuals:
##        1        2        3        4        5 
## -0.01250  0.01875  0.01875 -0.04375  0.01875 
## 
## Coefficients:
##                            Estimate Std. Error t value Pr(>|t|)    
## (Intercept)                 0.60625    0.01768  34.295 0.000849 ***
## poly(as.numeric(Sand), 2)1  0.09882    0.03953   2.500 0.129612    
## poly(as.numeric(Sand), 2)2 -0.23385    0.03953  -5.916 0.027402 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.03953 on 2 degrees of freedom
## Multiple R-squared:  0.9538, Adjusted R-squared:  0.9075 
## F-statistic: 20.62 on 2 and 2 DF,  p-value: 0.04624

Test of limiting factors - brome

recruit[is.na(recruit)] <- 0

## water
water <- subset(recruit, Treatment == "control" | Treatment=="water")

## germination
water.germ <- water %>% group_by(Density, Lvl) %>% summarize(eph=mean(ephedra.end ),brome=mean(brome.begin))

ggplot(water.germ) + geom_jitter(aes(x=Density, y=brome, fill=Lvl, color=Lvl), size=2)+ylab("number of Bromus germinants")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

m1 <- glm(brome.begin~Density *Lvl, data=water, family=poisson)
anova(m1, test="Chisq")
## Analysis of Deviance Table
## 
## Model: poisson, link: log
## 
## Response: brome.begin
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)    
## NULL                          299    1275.05             
## Density      1   849.52       298     425.53   <2e-16 ***
## Lvl          2     4.50       296     421.03   0.1053    
## Density:Lvl  2     0.00       294     421.03   0.9985    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## biomass
water.bio <- water %>% group_by(Density, Lvl) %>% summarize(eph=mean(Ephedra.biomass ),brome=mean(brome.biomass))

ggplot(water.bio) + geom_jitter(aes(x=Density, y=brome, fill=Lvl, color=Lvl), size=2)+ylab("final biomass of Brome")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

## clipping
clip <- subset(recruit, Treatment == "control" | Treatment=="clipped")

## germination
clip.germ <- clip %>% group_by(Density, Lvl) %>% summarize(eph=mean(ephedra.end ),brome=mean(brome.begin))

ggplot(clip.germ) + geom_jitter(aes(x=Density, y=brome, fill=Lvl, color=Lvl), size=2)+ylab("number of Bromus germinants")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

## biomass
clip.bio <- clip %>% group_by(Density, Lvl) %>% summarize(eph=mean(Ephedra.biomass ),brome=mean(brome.biomass))

ggplot(clip.bio) + geom_jitter(aes(x=Density, y=brome, fill=Lvl, color=Lvl), size=2)+ylab("final biomass of Brome")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

m1 <- glm(brome.begin~Density *Lvl, data=clip, family=poisson)
anova(m1, test="Chisq")
## Analysis of Deviance Table
## 
## Model: poisson, link: log
## 
## Response: brome.begin
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)    
## NULL                          299    1332.93             
## Density      1   863.30       298     469.63   <2e-16 ***
## Lvl          2     3.66       296     465.97   0.1605    
## Density:Lvl  2     0.04       294     465.93   0.9814    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## shade
shade <- subset(recruit, Treatment == "control" | Treatment=="shade")

## germination
shade.germ <- shade %>% group_by(Density, Lvl) %>% summarize(eph=mean(ephedra.end ),brome=mean(brome.begin))

ggplot(shade.germ) + geom_jitter(aes(x=Density, y=brome, fill=Lvl, color=Lvl), size=2)+ylab("number of Bromus germinants")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

## biomass
shade.bio <- shade %>% group_by(Density, Lvl) %>% summarize(eph=mean(Ephedra.biomass ),brome=mean(brome.biomass))

ggplot(shade.bio) + geom_jitter(aes(x=Density, y=brome, fill=Lvl, color=Lvl), size=2)+ylab("final biomass of Brome")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

m1 <- glm(brome.begin~Density *Lvl, data=shade, family=poisson)
anova(m1, test="Chisq")
## Analysis of Deviance Table
## 
## Model: poisson, link: log
## 
## Response: brome.begin
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
## NULL                          299    1221.52              
## Density      1   804.97       298     416.55 < 2.2e-16 ***
## Lvl          2    13.34       296     403.21  0.001267 ** 
## Density:Lvl  2     0.22       294     402.98  0.894507    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Test of limiting factors - Ephedra

## water
## germination

ggplot(water.germ) + geom_jitter(aes(x=Density, y=eph, fill=Lvl, color=Lvl), size=2)+ylab("number of Ephedra germinants")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

## biomass

ggplot(water.bio) + geom_jitter(aes(x=Density, y=eph, fill=Lvl, color=Lvl), size=2)+ylab("final biomass of Ephedra")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

## clipping

## germination

ggplot(clip.germ) + geom_jitter(aes(x=Density, y=eph, fill=Lvl, color=Lvl), size=2)+ylab("number of Ephedra germinants")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

## biomass

ggplot(clip.bio) + geom_jitter(aes(x=Density, y=eph, fill=Lvl, color=Lvl), size=2)+ylab("final biomass of Ephedra")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

m1 <- glm(ephedra.end~Density *Lvl, data=clip, family=poisson)
anova(m1, test="Chisq")
## Analysis of Deviance Table
## 
## Model: poisson, link: log
## 
## Response: ephedra.end
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)  
## NULL                          299     220.71           
## Density      1  2.73641       298     217.97  0.09808 .
## Lvl          2  0.34733       296     217.63  0.84058  
## Density:Lvl  2  0.47792       294     217.15  0.78745  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## shade

## germination

ggplot(shade.germ) + geom_jitter(aes(x=Density, y=eph, fill=Lvl, color=Lvl), size=2)+ylab("number of Ephedra germinants")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

##biomass
ggplot(shade.bio) + geom_jitter(aes(x=Density, y=eph, fill=Lvl, color=Lvl), size=2)+ylab("final biomass of Ephedra")+ theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),axis.line = element_line(colour = "black"),panel.background = element_blank(),axis.text=element_text(size=14),axis.title=element_text(size=16))

m1 <- glm(brome.begin~Density *Lvl, data=shade, family=poisson)
anova(m1, test="Chisq")
## Analysis of Deviance Table
## 
## Model: poisson, link: log
## 
## Response: brome.begin
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
## NULL                          299    1221.52              
## Density      1   804.97       298     416.55 < 2.2e-16 ***
## Lvl          2    13.34       296     403.21  0.001267 ** 
## Density:Lvl  2     0.22       294     402.98  0.894507    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Summarized GLMs for manuscript

## Water
m1 <- glm(ephedra.end~ Density *Lvl, data=water, family=binomial)
anova(m1, test="Chisq")
## Analysis of Deviance Table
## 
## Model: binomial, link: logit
## 
## Response: ephedra.end
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)  
## NULL                          299     397.45           
## Density      1  2.78375       298     394.66  0.09522 .
## Lvl          2  0.20096       296     394.46  0.90440  
## Density:Lvl  2  0.08676       294     394.37  0.95755  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
m1.bio <- glm(Ephedra.biomass~Density *Lvl, data=subset(water, Ephedra.biomass>0 & Ephedra.biomass<0.6), family=Gamma)
anova(m1.bio, test="Chisq")
## Analysis of Deviance Table
## 
## Model: Gamma, link: inverse
## 
## Response: Ephedra.biomass
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)
## NULL                          106     39.539         
## Density      1  0.11078       105     39.428   0.5364
## Lvl          2  0.10032       103     39.328   0.8411
## Density:Lvl  2  0.53107       101     38.796   0.4001
## Shade
m2 <- glm(ephedra.end~ Density *Lvl, data=shade, family=binomial)
anova(m2, test="Chisq")
## Analysis of Deviance Table
## 
## Model: binomial, link: logit
## 
## Response: ephedra.end
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)   
## NULL                          299     368.20            
## Density      1   0.7330       298     367.46 0.391927   
## Lvl          2  11.4709       296     355.99 0.003229 **
## Density:Lvl  2   1.4327       294     354.56 0.488533   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
m2.bio <- glm(Ephedra.biomass~Density *Lvl, data=subset(shade, Ephedra.biomass>0 & Ephedra.biomass < 0.7), family=Gamma)
anova(m2.bio, test="Chisq")
## Analysis of Deviance Table
## 
## Model: Gamma, link: inverse
## 
## Response: Ephedra.biomass
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)  
## NULL                           67     25.637           
## Density      1  0.39889        66     25.238  0.23231  
## Lvl          2  2.02160        64     23.217  0.02691 *
## Density:Lvl  2  0.02846        62     23.188  0.95038  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Clipping
m3 <- glm(ephedra.end~ Density *Lvl, data=clip, family=binomial)
anova(m3, test="Chisq")
## Analysis of Deviance Table
## 
## Model: binomial, link: logit
## 
## Response: ephedra.end
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)  
## NULL                          299     393.19           
## Density      1   4.2054       298     388.98   0.0403 *
## Lvl          2   0.5545       296     388.43   0.7579  
## Density:Lvl  2   0.8190       294     387.61   0.6640  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
m3.bio <- glm(Ephedra.biomass~Density *Lvl, data=subset(clip, Ephedra.biomass>0 & Ephedra.biomass < 0.7), family=Gamma)
anova(m3.bio, test="Chisq")
## Analysis of Deviance Table
## 
## Model: Gamma, link: inverse
## 
## Response: Ephedra.biomass
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)    
## NULL                          106     48.613             
## Density      1   0.6174       105     47.996   0.1774    
## Lvl          2   9.2947       103     38.701 1.13e-06 ***
## Density:Lvl  2   0.9905       101     37.711   0.2324    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
m4 <-  glm(ephedra.end~brome.begin, data=recruit, family=binomial)
anova(m4, test="Chisq")
## Analysis of Deviance Table
## 
## Model: binomial, link: logit
## 
## Response: ephedra.end
## 
## Terms added sequentially (first to last)
## 
## 
##             Df Deviance Resid. Df Resid. Dev Pr(>Chi)
## NULL                          699     901.37         
## brome.begin  1   2.0214       698     899.35   0.1551
m4.bio <-  glm(Ephedra.biomass~brome.begin, data=subset(recruit, Ephedra.biomass>0 & Ephedra.biomass < 0.7), family=Gamma)
anova(m4.bio, test="Chisq")
## Analysis of Deviance Table
## 
## Model: Gamma, link: inverse
## 
## Response: Ephedra.biomass
## 
## Terms added sequentially (first to last)
## 
## 
##             Df  Deviance Resid. Df Resid. Dev Pr(>Chi)
## NULL                           209     89.330         
## brome.begin  1 0.0012905       208     89.328    0.953
## brome comparisons
library(lme4)
## Loading required package: Matrix
## 
## Attaching package: 'Matrix'
## The following object is masked from 'package:tidyr':
## 
##     expand
## Water
m1 <- glm(brome.begin/Density~ Lvl, weights=Density, data=subset(water, Density>0), family=binomial)
anova(m1, test="Chisq")
## Analysis of Deviance Table
## 
## Model: binomial, link: logit
## 
## Response: brome.begin/Density
## 
## Terms added sequentially (first to last)
## 
## 
##      Df Deviance Resid. Df Resid. Dev Pr(>Chi)   
## NULL                   239     475.39            
## Lvl   2   10.831       237     464.56 0.004447 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Shade
m2 <- glm(brome.begin/Density~ Lvl, weights=Density, data=subset(shade, Density>0), family=binomial)
anova(m2, test="Chisq")
## Analysis of Deviance Table
## 
## Model: binomial, link: logit
## 
## Response: brome.begin/Density
## 
## Terms added sequentially (first to last)
## 
## 
##      Df Deviance Resid. Df Resid. Dev  Pr(>Chi)    
## NULL                   239     429.07              
## Lvl   2   31.496       237     397.57 1.448e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Clipping
m3<- glm(brome.begin/Density~ Lvl, weights=Density, data=subset(clip, Density>0), family=binomial)
anova(m3, test="Chisq")
## Analysis of Deviance Table
## 
## Model: binomial, link: logit
## 
## Response: brome.begin/Density
## 
## Terms added sequentially (first to last)
## 
## 
##      Df Deviance Resid. Df Resid. Dev Pr(>Chi)  
## NULL                   239     528.32           
## Lvl   2    8.962       237     519.35  0.01132 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Brome density vs Ephedra

dense <- recruit %>% filter(Ephedra.biomass>0) %>%  group_by(brome.begin) %>% summarize(bio=mean(Ephedra.biomass), recruit=mean(ephedra.end), germ=mean(ephedra.begin), abv.grd=mean(Ephedra.above),blw.grd=mean(Ephedra.below))

ggplot(dense, aes(x=brome.begin, y=bio))+geom_point()+ theme_Publication() +ylab("E. californica biomass") + xlab("brome density") + geom_smooth(method=lm, formula= y ~ poly(x,2))
## Warning: `legend.margin` must be specified using `margin()`. For the old
## behavior use legend.spacing
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x
## $y, : font family not found in Windows font database

m1 <- lm(bio~poly(brome.begin,2), data=dense)

## recruitment no effect
ggplot(dense, aes(x=brome.begin, y=recruit))+geom_point()+ theme_Publication() +ylab("E. californica biomass") + xlab("brome density") #+ geom_smooth(method=lm, formula= y ~ x)
## Warning: `legend.margin` must be specified using `margin()`. For the old
## behavior use legend.spacing
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x
## $y, : font family not found in Windows font database

m2 <- lm(recruit~brome.begin, data=dense)

## germination no effect
ggplot(dense, aes(x=brome.begin, y=germ))+geom_point()+ theme_Publication() +ylab("E. californica biomass") + xlab("brome density") #+ geom_smooth(method=lm, formula= y ~ x)
## Warning: `legend.margin` must be specified using `margin()`. For the old
## behavior use legend.spacing
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x
## $y, : font family not found in Windows font database

m3 <- lm(germ~brome.begin, data=dense)


## above ground biomass 
ggplot(dense, aes(x=brome.begin, y=abv.grd))+geom_point()+ theme_Publication() +ylab("above-ground biomass") + xlab("B. madritensis") + geom_smooth(method=lm, formula= y ~ x)
## Warning: `legend.margin` must be specified using `margin()`. For the old
## behavior use legend.spacing
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x
## $y, : font family not found in Windows font database

m4 <- lm(abv.grd~brome.begin, data=dense)
summary(m4)
## 
## Call:
## lm(formula = abv.grd ~ brome.begin, data = dense)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.010818 -0.003872 -0.001258  0.004408  0.013263 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  0.0365715  0.0034310  10.659 4.21e-08 ***
## brome.begin -0.0009681  0.0003850  -2.515   0.0247 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007263 on 14 degrees of freedom
## Multiple R-squared:  0.3112, Adjusted R-squared:  0.262 
## F-statistic: 6.325 on 1 and 14 DF,  p-value: 0.02474
## below ground biomass
ggplot(dense, aes(x=brome.begin, y=blw.grd))+geom_point()+ theme_Publication() +ylab("below-ground biomass") + xlab("B. madritensis density") + geom_smooth(method=lm, formula= y ~ poly(x,2))
## Warning: `legend.margin` must be specified using `margin()`. For the old
## behavior use legend.spacing
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database

## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, :
## font family not found in Windows font database
## Warning in grid.Call.graphics(C_text, as.graphicsAnnot(x$label), x$x, x
## $y, : font family not found in Windows font database

m5 <- lm(blw.grd~poly(brome.begin,2), data=dense)
summary(m5)
## 
## Call:
## lm(formula = blw.grd ~ poly(brome.begin, 2), data = dense)
## 
## Residuals:
##        Min         1Q     Median         3Q        Max 
## -0.0014710 -0.0007298 -0.0001058  0.0005620  0.0015789 
## 
## Coefficients:
##                         Estimate Std. Error t value Pr(>|t|)    
## (Intercept)            0.0049697  0.0002432   20.44 2.89e-11 ***
## poly(brome.begin, 2)1 -0.0036473  0.0009726   -3.75  0.00243 ** 
## poly(brome.begin, 2)2 -0.0032968  0.0009726   -3.39  0.00484 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.0009726 on 13 degrees of freedom
## Multiple R-squared:  0.6628, Adjusted R-squared:  0.6109 
## F-statistic: 12.78 on 2 and 13 DF,  p-value: 0.0008538

Germination patterns

recruit
##      ID Treatment   Lvl Rep Density ephedra.emergence branching fluoro
## 1     5   clipped  once   1       0                 0         0 0.0000
## 2    10   clipped  once   2       0                 1         2 0.3215
## 3    13   clipped  once   3       0                 1         6 0.7433
## 4    17   clipped  once   4       0                 0         0 0.0000
## 5    23   clipped  once   5       0                 0         0 0.0000
## 6    26   clipped  once   6       0                 0         0 0.0000
## 7    35   clipped  once   7       0                 0         0 0.0000
## 8    39   clipped  once   8       0                 1         4 0.6366
## 9    44   clipped  once   9       0                 0         0 0.0000
## 10   48   clipped  once  10       0                 0         0 0.0000
## 11   53   clipped  once  11       0                 2         8 0.7806
## 12   56   clipped  once  12       0                 0         0 0.0000
## 13   64   clipped  once  13       0                 0         0 0.0000
## 14   69   clipped  once  14       0                 0         0 0.0000
## 15   74   clipped  once  15       0                 0         0 0.0000
## 16   76   clipped  once  16       0                 0         0 0.0000
## 17   81   clipped  once  17       0                 0         0 0.0000
## 18   86   clipped  once  18       0                 0         0 0.0000
## 19   92   clipped  once  19       0                 0         0 0.0000
## 20   96   clipped  once  20       0                 0         0 0.0000
## 21  105   clipped twice   1       0                 0         0 0.0000
## 22  110   clipped twice   2       0                 1         8 0.7619
## 23  113   clipped twice   3       0                 1         5 0.8012
## 24  117   clipped twice   4       0                 1         4 0.7447
## 25  123   clipped twice   5       0                 2         5 0.7471
## 26  126   clipped twice   6       0                 0         0 0.0000
## 27  135   clipped twice   7       0                 1         4 0.5985
## 28  139   clipped twice   8       0                 1         5 0.7935
## 29  144   clipped twice   9       0                 0         0 0.0000
## 30  148   clipped twice  10       0                 1         6 0.6318
## 31  153   clipped twice  11       0                 1         6 0.6807
## 32  156   clipped twice  12       0                 0         0 0.0000
## 33  164   clipped twice  13       0                 0         0 0.0000
## 34  169   clipped twice  14       0                 0         0 0.0000
## 35  174   clipped twice  15       0                 0         0 0.0000
## 36  176   clipped twice  16       0                 1         4 0.7650
## 37  181   clipped twice  17       0                 0         0 0.0000
## 38  186   clipped twice  18       0                 0         0 0.0000
## 39  192   clipped twice  19       0                 0         0 0.0000
## 40  196   clipped twice  20       0                 0         0 0.0000
## 41    4   clipped  once   1       2                 0         0 0.0000
## 42    9   clipped  once   2       2                 0         0 0.0000
## 43   15   clipped  once   3       2                 0         0 0.0000
## 44   20   clipped  once   4       2                 0         0 0.0000
## 45   21   clipped  once   5       2                 0         0 0.0000
## 46   28   clipped  once   6       2                 0         0 0.0000
## 47   34   clipped  once   7       2                 0         0 0.0000
## 48   38   clipped  once   8       2                 0         0 0.0000
## 49   42   clipped  once   9       2                 0         0 0.0000
## 50   47   clipped  once  10       2                 0         0 0.0000
## 51   52   clipped  once  11       2                 0         0 0.0000
## 52   60   clipped  once  12       2                 0         0 0.0000
## 53   63   clipped  once  13       2                 0         0 0.0000
## 54   70   clipped  once  14       2                 0         0 0.0000
## 55   71   clipped  once  15       2                 1         6 0.7615
## 56   80   clipped  once  16       2                 2         7 0.7943
## 57   83   clipped  once  17       2                 0         0 0.0000
## 58   88   clipped  once  18       2                 0         0 0.0000
## 59   91   clipped  once  19       2                 0         0 0.0000
## 60   98   clipped  once  20       2                 0         0 0.0000
## 61  104   clipped twice   1       2                 1         6 0.6954
## 62  109   clipped twice   2       2                 0         0 0.0000
## 63  115   clipped twice   3       2                 1         4 0.7570
## 64  120   clipped twice   4       2                 0         0 0.0000
## 65  121   clipped twice   5       2                 0         0 0.0000
## 66  128   clipped twice   6       2                 0         0 0.0000
## 67  134   clipped twice   7       2                 0         0 0.0000
## 68  138   clipped twice   8       2                 0         0 0.0000
## 69  142   clipped twice   9       2                 0         0 0.0000
## 70  147   clipped twice  10       2                 0         0 0.0000
## 71  152   clipped twice  11       2                 1         4 0.7681
## 72  160   clipped twice  12       2                 0         0 0.0000
## 73  163   clipped twice  13       2                 0         0 0.0000
## 74  170   clipped twice  14       2                 1         6 0.7063
## 75  171   clipped twice  15       2                 0         0 0.0000
## 76  180   clipped twice  16       2                 1         4 0.7995
## 77  183   clipped twice  17       2                 0         0 0.0000
## 78  188   clipped twice  18       2                 0         0 0.0000
## 79  191   clipped twice  19       2                 0         0 0.0000
## 80  198   clipped twice  20       2                 0         0 0.0000
## 81    3   clipped  once   1       5                 1         3 0.6416
## 82    8   clipped  once   2       5                 0         0 0.0000
## 83   14   clipped  once   3       5                 1         4 0.7920
## 84   18   clipped  once   4       5                 0         0 0.0000
## 85   22   clipped  once   5       5                 0         0 0.0000
## 86   29   clipped  once   6       5                 1         5 0.6169
## 87   32   clipped  once   7       5                 0         0 0.0000
## 88   40   clipped  once   8       5                 1         4 0.6732
## 89   41   clipped  once   9       5                 2         8 0.7775
## 90   49   clipped  once  10       5                 1         5 0.8054
## 91   54   clipped  once  11       5                 2         4 0.7938
## 92   59   clipped  once  12       5                 1         3 0.7967
## 93   62   clipped  once  13       5                 1         5 0.7612
## 94   66   clipped  once  14       5                 1         6 0.7740
## 95   73   clipped  once  15       5                 1        11 0.7853
## 96   77   clipped  once  16       5                 0         0 0.0000
## 97   85   clipped  once  17       5                 0         0 0.0000
## 98   87   clipped  once  18       5                 1         6 0.7532
## 99   95   clipped  once  19       5                 0         0 0.0000
## 100  97   clipped  once  20       5                 1         5 0.7839
## 101 103   clipped twice   1       5                 2         6 0.7245
## 102 108   clipped twice   2       5                 2         6 0.7800
## 103 114   clipped twice   3       5                 0         0 0.0000
## 104 118   clipped twice   4       5                 0         0 0.0000
## 105 122   clipped twice   5       5                 1         4 0.7698
## 106 129   clipped twice   6       5                 0         0 0.0000
## 107 132   clipped twice   7       5                 1         6 0.7884
## 108 140   clipped twice   8       5                 0         0 0.0000
## 109 141   clipped twice   9       5                 0         0 0.0000
## 110 149   clipped twice  10       5                 0         0 0.0000
## 111 154   clipped twice  11       5                 0         0 0.0000
## 112 159   clipped twice  12       5                 0         0 0.0000
## 113 162   clipped twice  13       5                 1         4 0.4452
## 114 166   clipped twice  14       5                 1         7 0.6691
## 115 173   clipped twice  15       5                 0         0 0.0000
## 116 177   clipped twice  16       5                 0         0 0.0000
## 117 185   clipped twice  17       5                 1         4 0.7774
## 118 187   clipped twice  18       5                 0         0 0.0000
## 119 195   clipped twice  19       5                 1         5 0.7694
## 120 197   clipped twice  20       5                 1         4 0.7631
## 121   1   clipped  once   1      10                 1         6 0.7867
## 122   6   clipped  once   2      10                 1         3 0.7592
## 123  12   clipped  once   3      10                 0         0 0.0000
## 124  16   clipped  once   4      10                 0         0 0.0000
## 125  25   clipped  once   5      10                 0         0 0.0000
## 126  27   clipped  once   6      10                 0         0 0.0000
## 127  33   clipped  once   7      10                 1         5 0.7789
## 128  36   clipped  once   8      10                 1         7 0.7092
## 129  43   clipped  once   9      10                 0         0 0.0000
## 130  50   clipped  once  10      10                 0         0 0.0000
## 131  55   clipped  once  11      10                 1         4 0.6851
## 132  58   clipped  once  12      10                 0         0 0.0000
## 133  65   clipped  once  13      10                 0         0 0.0000
## 134  68   clipped  once  14      10                 0         0 0.0000
## 135  75   clipped  once  15      10                 1         4 0.7576
## 136  78   clipped  once  16      10                 1         8 0.7499
## 137  84   clipped  once  17      10                 0         0 0.0000
## 138  90   clipped  once  18      10                 0         0 0.0000
## 139  93   clipped  once  19      10                 0         0 0.0000
## 140  99   clipped  once  20      10                 0         0 0.0000
## 141 101   clipped twice   1      10                 0         0 0.0000
## 142 106   clipped twice   2      10                 1         6 0.7842
## 143 112   clipped twice   3      10                 0         0 0.0000
## 144 116   clipped twice   4      10                 0         0 0.0000
## 145 125   clipped twice   5      10                 0         0 0.0000
## 146 127   clipped twice   6      10                 1         5 0.7753
## 147 133   clipped twice   7      10                 0         0 0.0000
## 148 136   clipped twice   8      10                 0         0 0.0000
## 149 143   clipped twice   9      10                 0         0 0.0000
## 150 150   clipped twice  10      10                 1         5 0.7420
## 151 155   clipped twice  11      10                 0         0 0.0000
## 152 158   clipped twice  12      10                 2         7 0.7999
## 153 165   clipped twice  13      10                 0         0 0.0000
## 154 168   clipped twice  14      10                 3         9 0.6319
## 155 175   clipped twice  15      10                 0         0 0.0000
## 156 178   clipped twice  16      10                 2         5 0.5698
## 157 184   clipped twice  17      10                 0         0 0.0000
## 158 190   clipped twice  18      10                 1         4 0.7780
## 159 193   clipped twice  19      10                 0         0 0.0000
## 160 199   clipped twice  20      10                 0         0 0.0000
## 161   2   clipped  once   1      20                 1         7 0.7863
## 162   7   clipped  once   2      20                 0         0 0.0000
## 163  11   clipped  once   3      20                 0         0 0.0000
## 164  19   clipped  once   4      20                 0         0 0.0000
## 165  24   clipped  once   5      20                 0         0 0.0000
## 166  30   clipped  once   6      20                 0         0 0.0000
## 167  31   clipped  once   7      20                 0         0 0.0000
## 168  37   clipped  once   8      20                 0         0 0.0000
## 169  45   clipped  once   9      20                 0         0 0.0000
## 170  46   clipped  once  10      20                 0         0 0.0000
## 171  51   clipped  once  11      20                 0         0 0.0000
## 172  57   clipped  once  12      20                 0         0 0.0000
## 173  61   clipped  once  13      20                 0         0 0.0000
## 174  67   clipped  once  14      20                 0         0 0.0000
## 175  72   clipped  once  15      20                 0         0 0.0000
## 176  79   clipped  once  16      20                 0         0 0.0000
## 177  82   clipped  once  17      20                 0         0 0.0000
## 178  89   clipped  once  18      20                 0         0 0.0000
## 179  94   clipped  once  19      20                 0         0 0.0000
## 180 100   clipped  once  20      20                 1         4 0.6715
## 181 102   clipped twice   1      20                 0         0 0.0000
## 182 107   clipped twice   2      20                 0         0 0.0000
## 183 111   clipped twice   3      20                 0         0 0.0000
## 184 119   clipped twice   4      20                 0         0 0.0000
## 185 124   clipped twice   5      20                 0         0 0.0000
## 186 130   clipped twice   6      20                 1         5 0.5799
## 187 131   clipped twice   7      20                 0         0 0.0000
## 188 137   clipped twice   8      20                 0         0 0.0000
## 189 145   clipped twice   9      20                 1         3 0.7009
## 190 146   clipped twice  10      20                 0         0 0.0000
## 191 151   clipped twice  11      20                 1         5 0.7652
## 192 157   clipped twice  12      20                 0         0 0.0000
## 193 161   clipped twice  13      20                 0         0 0.0000
## 194 167   clipped twice  14      20                 0         0 0.0000
## 195 172   clipped twice  15      20                 0         0 0.0000
## 196 179   clipped twice  16      20                 0         0 0.0000
## 197 182   clipped twice  17      20                 2         6 0.7488
## 198 189   clipped twice  18      20                 0         0 0.0000
## 199 194   clipped twice  19      20                 0         0 0.0000
## 200 200   clipped twice  20      20                 0         0 0.0000
## 201 205   control  none   1       0                 0         0 0.0000
## 202 210   control  none   2       0                 0         0 0.0000
## 203 213   control  none   3       0                 1         5 0.7908
## 204 217   control  none   4       0                 1         3 0.7910
## 205 223   control  none   5       0                 0         0 0.0000
## 206 226   control  none   6       0                 0         0 0.0000
## 207 235   control  none   7       0                 1         6 0.7779
## 208 239   control  none   8       0                 0         0 0.0000
## 209 244   control  none   9       0                 1         6 0.5830
## 210 248   control  none  10       0                 1         5 0.7768
## 211 253   control  none  11       0                 0         0 0.0000
## 212 256   control  none  12       0                 0         0 0.0000
## 213 264   control  none  13       0                 0         0 0.0000
## 214 269   control  none  14       0                 0         0 0.0000
## 215 274   control  none  15       0                 0         0 0.0000
## 216 276   control  none  16       0                 1         4 0.7781
## 217 281   control  none  17       0                 0         0 0.0000
## 218 286   control  none  18       0                 0         0 0.0000
## 219 292   control  none  19       0                 0         0 0.0000
## 220 296   control  none  20       0                 0         0 0.0000
## 221 204   control  none   1       2                 1         9 0.6990
## 222 209   control  none   2       2                 0         0 0.0000
## 223 215   control  none   3       2                 0         0 0.0000
## 224 220   control  none   4       2                 0         0 0.0000
## 225 221   control  none   5       2                 0         0 0.0000
## 226 228   control  none   6       2                 0         0 0.0000
## 227 234   control  none   7       2                 1         5 0.7909
## 228 238   control  none   8       2                 0         0 0.0000
## 229 242   control  none   9       2                 1         5 0.7632
## 230 247   control  none  10       2                 0         0 0.0000
## 231 252   control  none  11       2                 1         7 0.8060
## 232 260   control  none  12       2                 1         7 0.7774
## 233 263   control  none  13       2                 0         0 0.0000
## 234 270   control  none  14       2                 0         0 0.0000
## 235 271   control  none  15       2                 1         6 0.7727
## 236 280   control  none  16       2                 0         0 0.0000
## 237 283   control  none  17       2                 2         5 0.7937
## 238 288   control  none  18       2                 0         0 0.0000
## 239 291   control  none  19       2                 1         3 0.4486
## 240 298   control  none  20       2                 0         0 0.0000
## 241 203   control  none   1       5                 0         0 0.0000
## 242 208   control  none   2       5                 0         0 0.0000
## 243 214   control  none   3       5                 0         0 0.0000
## 244 218   control  none   4       5                 0         0 0.0000
## 245 222   control  none   5       5                 0         0 0.0000
## 246 229   control  none   6       5                 1         6 0.8044
## 247 232   control  none   7       5                 1         3 0.8003
## 248 240   control  none   8       5                 0         0 0.0000
## 249 241   control  none   9       5                 1         3 0.4088
## 250 249   control  none  10       5                 0         0 0.0000
## 251 254   control  none  11       5                 0         0 0.0000
## 252 259   control  none  12       5                 0         0 0.0000
## 253 262   control  none  13       5                 1         9 0.7871
## 254 266   control  none  14       5                 2         5 0.7579
## 255 273   control  none  15       5                 0         0 0.0000
## 256 277   control  none  16       5                 1         3 0.7850
## 257 285   control  none  17       5                 2         8 0.7973
## 258 287   control  none  18       5                 0         0 0.0000
## 259 295   control  none  19       5                 1         5 0.7066
## 260 297   control  none  20       5                 0         0 0.0000
## 261 201   control  none   1      10                 0         0 0.0000
## 262 206   control  none   2      10                 1         5 0.7773
## 263 212   control  none   3      10                 0         0 0.0000
## 264 216   control  none   4      10                 0         0 0.0000
## 265 225   control  none   5      10                 0         0 0.0000
## 266 227   control  none   6      10                 1         7 0.7974
## 267 233   control  none   7      10                 2         3 0.6053
## 268 236   control  none   8      10                 0         0 0.0000
## 269 243   control  none   9      10                 0         0 0.0000
## 270 250   control  none  10      10                 0         0 0.0000
## 271 255   control  none  11      10                 2         7 0.7887
## 272 258   control  none  12      10                 0         0 0.0000
## 273 265   control  none  13      10                 0         0 0.0000
## 274 268   control  none  14      10                 0         0 0.0000
## 275 275   control  none  15      10                 1         7 0.7880
## 276 278   control  none  16      10                 0         0 0.0000
## 277 284   control  none  17      10                 1         4 0.7938
## 278 290   control  none  18      10                 0         0 0.0000
## 279 293   control  none  19      10                 0         0 0.0000
## 280 299   control  none  20      10                 0         0 0.0000
## 281 202   control  none   1      20                 1         3 0.7360
## 282 207   control  none   2      20                 1         5 0.7544
## 283 211   control  none   3      20                 0         0 0.0000
## 284 219   control  none   4      20                 0         0 0.0000
## 285 224   control  none   5      20                 0         0 0.0000
## 286 230   control  none   6      20                 0         0 0.0000
## 287 231   control  none   7      20                 2         5 0.6865
## 288 237   control  none   8      20                 1         3 0.7434
## 289 245   control  none   9      20                 0         0 0.0000
## 290 246   control  none  10      20                 0         0 0.0000
## 291 251   control  none  11      20                 0         0 0.0000
## 292 257   control  none  12      20                 2         5 0.5452
## 293 261   control  none  13      20                 0         0 0.0000
## 294 267   control  none  14      20                 0         0 0.0000
## 295 272   control  none  15      20                 0         0 0.0000
## 296 279   control  none  16      20                 0         0 0.0000
## 297 282   control  none  17      20                 0         0 0.0000
## 298 289   control  none  18      20                 0         0 0.0000
## 299 294   control  none  19      20                 0         0 0.0000
## 300 300   control  none  20      20                 2         7 0.7728
## 301 305     shade  high   1       0                 0         0 0.0000
## 302 310     shade  high   2       0                 0         0 0.0000
## 303 313     shade  high   3       0                 0         0 0.0000
## 304 317     shade  high   4       0                 1         3 0.6200
## 305 323     shade  high   5       0                 0         0 0.0000
## 306 326     shade  high   6       0                 0         0 0.0000
## 307 335     shade  high   7       0                 0         0 0.0000
## 308 339     shade  high   8       0                 0         0 0.0000
## 309 344     shade  high   9       0                 0         0 0.0000
## 310 348     shade  high  10       0                 0         0 0.0000
## 311 353     shade  high  11       0                 0         0 0.0000
## 312 356     shade  high  12       0                 0         0 0.0000
## 313 364     shade  high  13       0                 0         0 0.0000
## 314 369     shade  high  14       0                 0         0 0.0000
## 315 374     shade  high  15       0                 0         0 0.0000
## 316 376     shade  high  16       0                 0         0 0.0000
## 317 381     shade  high  17       0                 0         0 0.0000
## 318 386     shade  high  18       0                 0         0 0.0000
## 319 392     shade  high  19       0                 0         0 0.0000
## 320 396     shade  high  20       0                 1         3 0.8874
## 321 405     shade   med   1       0                 0         0 0.0000
## 322 410     shade   med   2       0                 0         0 0.0000
## 323 413     shade   med   3       0                 0         0 0.0000
## 324 417     shade   med   4       0                 0         0 0.0000
## 325 423     shade   med   5       0                 0         0 0.0000
## 326 426     shade   med   6       0                 1         8 0.7693
## 327 435     shade   med   7       0                 1         7 0.7306
## 328 439     shade   med   8       0                 0         0 0.0000
## 329 444     shade   med   9       0                 0         0 0.0000
## 330 448     shade   med  10       0                 0         0 0.0000
## 331 453     shade   med  11       0                 0         0 0.0000
## 332 456     shade   med  12       0                 0         0 0.0000
## 333 464     shade   med  13       0                 0         0 0.0000
## 334 469     shade   med  14       0                 0         0 0.0000
## 335 474     shade   med  15       0                 0         0 0.0000
## 336 476     shade   med  16       0                 0         0 0.0000
## 337 481     shade   med  17       0                 0         0 0.0000
## 338 486     shade   med  18       0                 1         7 0.7961
## 339 492     shade   med  19       0                 0         0 0.0000
## 340 496     shade   med  20       0                 0         0 0.0000
## 341 304     shade  high   1       2                 0         0 0.0000
## 342 309     shade  high   2       2                 0         0 0.0000
## 343 315     shade  high   3       2                 0         0 0.0000
## 344 320     shade  high   4       2                 1         4 0.7831
## 345 321     shade  high   5       2                 1         5 0.8014
## 346 328     shade  high   6       2                 0         0 0.0000
## 347 334     shade  high   7       2                 0         0 0.0000
## 348 338     shade  high   8       2                 0         0 0.0000
## 349 342     shade  high   9       2                 0         0 0.0000
## 350 347     shade  high  10       2                 0         0 0.0000
## 351 352     shade  high  11       2                 0         0 0.0000
## 352 360     shade  high  12       2                 0         0 0.0000
## 353 363     shade  high  13       2                 2         5 0.7799
## 354 370     shade  high  14       2                 0         0 0.0000
## 355 371     shade  high  15       2                 0         0 0.0000
## 356 380     shade  high  16       2                 0         0 0.0000
## 357 383     shade  high  17       2                 1         5 0.7620
## 358 388     shade  high  18       2                 0         0 0.0000
## 359 391     shade  high  19       2                 0         0 0.0000
## 360 398     shade  high  20       2                 0         0 0.0000
## 361 404     shade   med   1       2                 0         0 0.0000
## 362 409     shade   med   2       2                 0         0 0.0000
## 363 415     shade   med   3       2                 0         0 0.0000
## 364 420     shade   med   4       2                 0         0 0.0000
## 365 421     shade   med   5       2                 0         0 0.0000
## 366 428     shade   med   6       2                 1         4 0.7959
## 367 434     shade   med   7       2                 0         0 0.0000
## 368 438     shade   med   8       2                 0         0 0.0000
## 369 442     shade   med   9       2                 0         0 0.0000
## 370 447     shade   med  10       2                 0         0 0.0000
## 371 452     shade   med  11       2                 0         0 0.0000
## 372 460     shade   med  12       2                 0         0 0.0000
## 373 463     shade   med  13       2                 0         0 0.0000
## 374 470     shade   med  14       2                 0         0 0.0000
## 375 471     shade   med  15       2                 1         5 0.5438
## 376 480     shade   med  16       2                 0         0 0.0000
## 377 483     shade   med  17       2                 0         0 0.0000
## 378 488     shade   med  18       2                 2         8 0.7381
## 379 491     shade   med  19       2                 0         0 0.0000
## 380 498     shade   med  20       2                 0         0 0.0000
## 381 303     shade  high   1       5                 1         4 0.5240
## 382 308     shade  high   2       5                 1         4 0.7704
## 383 314     shade  high   3       5                 0         0 0.0000
## 384 318     shade  high   4       5                 0         0 0.0000
## 385 322     shade  high   5       5                 0         0 0.0000
## 386 329     shade  high   6       5                 0         0 0.0000
## 387 332     shade  high   7       5                 0         0 0.0000
## 388 340     shade  high   8       5                 0         0 0.0000
## 389 341     shade  high   9       5                 0         0 0.0000
## 390 349     shade  high  10       5                 0         0 0.0000
## 391 354     shade  high  11       5                 0         0 0.0000
## 392 359     shade  high  12       5                 0         0 0.0000
## 393 362     shade  high  13       5                 0         0 0.0000
## 394 366     shade  high  14       5                 0         0 0.0000
## 395 373     shade  high  15       5                 0         0 0.0000
## 396 377     shade  high  16       5                 0         0 0.0000
## 397 385     shade  high  17       5                 0         0 0.0000
## 398 387     shade  high  18       5                 0         0 0.0000
## 399 395     shade  high  19       5                 0         0 0.0000
## 400 397     shade  high  20       5                 0         0 0.0000
## 401 403     shade   med   1       5                 0         0 0.0000
## 402 408     shade   med   2       5                 2         6 0.5096
## 403 414     shade   med   3       5                 0         0 0.0000
## 404 418     shade   med   4       5                 0         0 0.0000
## 405 422     shade   med   5       5                 0         0 0.0000
## 406 429     shade   med   6       5                 0         0 0.0000
## 407 432     shade   med   7       5                 0         0 0.0000
## 408 440     shade   med   8       5                 1         5 0.8030
## 409 441     shade   med   9       5                 0         0 0.0000
## 410 449     shade   med  10       5                 0         0 0.0000
## 411 454     shade   med  11       5                 0         0 0.0000
## 412 459     shade   med  12       5                 0         0 0.0000
## 413 462     shade   med  13       5                 1         7 0.7493
## 414 466     shade   med  14       5                 0         0 0.0000
## 415 473     shade   med  15       5                 0         0 0.0000
## 416 477     shade   med  16       5                 0         0 0.0000
## 417 485     shade   med  17       5                 0         0 0.0000
## 418 487     shade   med  18       5                 0         0 0.0000
## 419 495     shade   med  19       5                 0         0 0.0000
## 420 497     shade   med  20       5                 0         0 0.0000
## 421 301     shade  high   1      10                 1         5 0.7900
## 422 306     shade  high   2      10                 0         0 0.0000
## 423 312     shade  high   3      10                 0         0 0.0000
## 424 316     shade  high   4      10                 0         0 0.0000
## 425 325     shade  high   5      10                 0         0 0.0000
## 426 327     shade  high   6      10                 0         0 0.0000
## 427 333     shade  high   7      10                 0         0 0.0000
## 428 336     shade  high   8      10                 0         0 0.0000
## 429 343     shade  high   9      10                 0         0 0.0000
## 430 350     shade  high  10      10                 0         0 0.0000
## 431 355     shade  high  11      10                 0         0 0.0000
## 432 358     shade  high  12      10                 0         0 0.0000
## 433 365     shade  high  13      10                 0         0 0.0000
## 434 368     shade  high  14      10                 0         0 0.0000
## 435 375     shade  high  15      10                 0         0 0.0000
## 436 378     shade  high  16      10                 0         0 0.0000
## 437 384     shade  high  17      10                 0         0 0.0000
## 438 390     shade  high  18      10                 0         0 0.0000
## 439 393     shade  high  19      10                 0         0 0.0000
## 440 399     shade  high  20      10                 0         0 0.0000
## 441 401     shade   med   1      10                 1         7 0.7954
## 442 406     shade   med   2      10                 0         0 0.0000
## 443 412     shade   med   3      10                 0         0 0.0000
## 444 416     shade   med   4      10                 0         0 0.0000
## 445 425     shade   med   5      10                 0         0 0.0000
## 446 427     shade   med   6      10                 0         0 0.0000
## 447 433     shade   med   7      10                 0         0 0.0000
## 448 436     shade   med   8      10                 0         0 0.0000
## 449 443     shade   med   9      10                 0         0 0.0000
## 450 450     shade   med  10      10                 0         0 0.0000
## 451 455     shade   med  11      10                 0         0 0.0000
## 452 458     shade   med  12      10                 0         0 0.0000
## 453 465     shade   med  13      10                 0         0 0.0000
## 454 468     shade   med  14      10                 0         0 0.0000
## 455 475     shade   med  15      10                 0         0 0.0000
## 456 478     shade   med  16      10                 0         0 0.0000
## 457 484     shade   med  17      10                 0         0 0.0000
## 458 490     shade   med  18      10                 0         0 0.0000
## 459 493     shade   med  19      10                 0         0 0.0000
## 460 499     shade   med  20      10                 0         0 0.0000
## 461 302     shade  high   1      20                 0         0 0.0000
## 462 307     shade  high   2      20                 0         0 0.0000
## 463 311     shade  high   3      20                 0         0 0.0000
## 464 319     shade  high   4      20                 0         0 0.0000
## 465 324     shade  high   5      20                 0         0 0.0000
## 466 330     shade  high   6      20                 0         0 0.0000
## 467 331     shade  high   7      20                 0         0 0.0000
## 468 337     shade  high   8      20                 0         0 0.0000
## 469 345     shade  high   9      20                 0         0 0.0000
## 470 346     shade  high  10      20                 1         7 0.7807
## 471 351     shade  high  11      20                 1         4 0.7201
## 472 357     shade  high  12      20                 0         0 0.0000
## 473 361     shade  high  13      20                 0         0 0.0000
## 474 367     shade  high  14      20                 0         0 0.0000
## 475 372     shade  high  15      20                 0         0 0.0000
## 476 379     shade  high  16      20                 0         0 0.0000
## 477 382     shade  high  17      20                 0         0 0.0000
## 478 389     shade  high  18      20                 0         0 0.0000
## 479 394     shade  high  19      20                 0         0 0.0000
## 480 400     shade  high  20      20                 0         0 0.0000
## 481 402     shade   med   1      20                 0         0 0.0000
## 482 407     shade   med   2      20                 0         0 0.0000
## 483 411     shade   med   3      20                 0         0 0.0000
## 484 419     shade   med   4      20                 0         0 0.0000
## 485 424     shade   med   5      20                 0         0 0.0000
## 486 430     shade   med   6      20                 0         0 0.0000
## 487 431     shade   med   7      20                 0         0 0.0000
## 488 437     shade   med   8      20                 0         0 0.0000
## 489 445     shade   med   9      20                 0         0 0.0000
## 490 446     shade   med  10      20                 0         0 0.0000
## 491 451     shade   med  11      20                 0         0 0.0000
## 492 457     shade   med  12      20                 0         0 0.0000
## 493 461     shade   med  13      20                 1         2 0.6052
## 494 467     shade   med  14      20                 0         0 0.0000
## 495 472     shade   med  15      20                 0         0 0.0000
## 496 479     shade   med  16      20                 0         0 0.0000
## 497 482     shade   med  17      20                 0         0 0.0000
## 498 489     shade   med  18      20                 0         0 0.0000
## 499 494     shade   med  19      20                 0         0 0.0000
## 500 500     shade   med  20      20                 0         0 0.0000
## 501 505     water  high   1       0                 1         5 0.7860
## 502 510     water  high   2       0                 1         5 0.7988
## 503 513     water  high   3       0                 1         4 0.8049
## 504 517     water  high   4       0                 2         7 0.7890
## 505 523     water  high   5       0                 0         0 0.0000
## 506 526     water  high   6       0                 1         6 0.7990
## 507 535     water  high   7       0                 0         0 0.0000
## 508 539     water  high   8       0                 0         0 0.0000
## 509 544     water  high   9       0                 1         5 0.8016
## 510 548     water  high  10       0                 0         0 0.0000
## 511 553     water  high  11       0                 0         0 0.0000
## 512 556     water  high  12       0                 0         0 0.0000
## 513 564     water  high  13       0                 1         4 0.3523
## 514 569     water  high  14       0                 0         0 0.0000
## 515 574     water  high  15       0                 1         7 0.7320
## 516 576     water  high  16       0                 0         0 0.0000
## 517 581     water  high  17       0                 0         0 0.0000
## 518 586     water  high  18       0                 0         0 0.0000
## 519 592     water  high  19       0                 1         6 0.8004
## 520 596     water  high  20       0                 0         0 0.0000
## 521 605     water   low   1       0                 1         7 0.8077
## 522 610     water   low   2       0                 0         0 0.0000
## 523 613     water   low   3       0                 1         5 0.7018
## 524 617     water   low   4       0                 1         8 0.7353
## 525 623     water   low   5       0                 0         0 0.0000
## 526 626     water   low   6       0                 1         8 0.7515
## 527 635     water   low   7       0                 1         7 0.7973
## 528 639     water   low   8       0                 0         0 0.0000
## 529 644     water   low   9       0                 0         0 0.0000
## 530 648     water   low  10       0                 1         5 0.7796
## 531 653     water   low  11       0                 0         0 0.0000
## 532 656     water   low  12       0                 0         0 0.0000
## 533 664     water   low  13       0                 0         0 0.0000
## 534 669     water   low  14       0                 2         9 0.7844
## 535 674     water   low  15       0                 0         0 0.0000
## 536 676     water   low  16       0                 0         0 0.0000
## 537 681     water   low  17       0                 0         0 0.0000
## 538 686     water   low  18       0                 0         0 0.0000
## 539 692     water   low  19       0                 0         0 0.0000
## 540 696     water   low  20       0                 0         0 0.0000
## 541 504     water  high   1       2                 0         0 0.0000
## 542 509     water  high   2       2                 0         0 0.0000
## 543 515     water  high   3       2                 0         0 0.0000
## 544 520     water  high   4       2                 0         0 0.0000
## 545 521     water  high   5       2                 1         6 0.5436
## 546 528     water  high   6       2                 0         0 0.0000
## 547 534     water  high   7       2                 1         4 0.6423
## 548 538     water  high   8       2                 0         0 0.0000
## 549 542     water  high   9       2                 1         6 0.7139
## 550 547     water  high  10       2                 0         0 0.0000
## 551 552     water  high  11       2                 1         3 0.6860
## 552 560     water  high  12       2                 0         0 0.0000
## 553 563     water  high  13       2                 0         0 0.0000
## 554 570     water  high  14       2                 0         0 0.0000
## 555 571     water  high  15       2                 0         0 0.0000
## 556 580     water  high  16       2                 2         4 0.5730
## 557 583     water  high  17       2                 1         6 0.7429
## 558 588     water  high  18       2                 1         6 0.7680
## 559 591     water  high  19       2                 1         5 0.6706
## 560 598     water  high  20       2                 0         0 0.0000
## 561 604     water   low   1       2                 1         6 0.6032
## 562 609     water   low   2       2                 0         0 0.0000
## 563 615     water   low   3       2                 0         0 0.0000
## 564 620     water   low   4       2                 0         0 0.0000
## 565 621     water   low   5       2                 0         0 0.0000
## 566 628     water   low   6       2                 1         3 0.7533
## 567 634     water   low   7       2                 0         0 0.0000
## 568 638     water   low   8       2                 0         0 0.0000
## 569 642     water   low   9       2                 1         3 0.1651
## 570 647     water   low  10       2                 1         9 0.7817
## 571 652     water   low  11       2                 0         0 0.0000
## 572 660     water   low  12       2                 0         0 0.0000
## 573 663     water   low  13       2                 1         7 0.7852
## 574 670     water   low  14       2                 1         6 0.7514
## 575 671     water   low  15       2                 0         0 0.0000
## 576 680     water   low  16       2                 0         0 0.0000
## 577 683     water   low  17       2                 1         5 0.7886
## 578 688     water   low  18       2                 0         0 0.0000
## 579 691     water   low  19       2                 0         0 0.0000
## 580 698     water   low  20       2                 1         5 0.8047
## 581 503     water  high   1       5                 1         5 0.7532
## 582 508     water  high   2       5                 0         0 0.0000
## 583 514     water  high   3       5                 1         9 0.7921
## 584 518     water  high   4       5                 1         7 0.7809
## 585 522     water  high   5       5                 1         4 0.7513
## 586 529     water  high   6       5                 0         0 0.0000
## 587 532     water  high   7       5                 0         0 0.0000
## 588 540     water  high   8       5                 0         0 0.0000
## 589 541     water  high   9       5                 0         0 0.0000
## 590 549     water  high  10       5                 0         0 0.0000
## 591 554     water  high  11       5                 0         0 0.0000
## 592 559     water  high  12       5                 0         0 0.0000
## 593 562     water  high  13       5                 0         0 0.0000
## 594 566     water  high  14       5                 0         0 0.0000
## 595 573     water  high  15       5                 0         0 0.0000
## 596 577     water  high  16       5                 0         0 0.0000
## 597 585     water  high  17       5                 0         0 0.0000
## 598 587     water  high  18       5                 0         0 0.0000
## 599 595     water  high  19       5                 0         0 0.0000
## 600 597     water  high  20       5                 0         0 0.0000
## 601 603     water   low   1       5                 0         0 0.0000
## 602 608     water   low   2       5                 0         0 0.0000
## 603 614     water   low   3       5                 0         0 0.0000
## 604 618     water   low   4       5                 0         0 0.0000
## 605 622     water   low   5       5                 0         0 0.0000
## 606 629     water   low   6       5                 1         5 0.7865
## 607 632     water   low   7       5                 0         0 0.0000
## 608 640     water   low   8       5                 0         0 0.0000
## 609 641     water   low   9       5                 1         5 0.7016
## 610 649     water   low  10       5                 0         0 0.0000
## 611 654     water   low  11       5                 0         0 0.0000
## 612 659     water   low  12       5                 0         0 0.0000
## 613 662     water   low  13       5                 0         0 0.0000
## 614 666     water   low  14       5                 0         0 0.0000
## 615 673     water   low  15       5                 0         0 0.0000
## 616 677     water   low  16       5                 0         0 0.0000
## 617 685     water   low  17       5                 1         5 0.7909
## 618 687     water   low  18       5                 0         0 0.0000
## 619 695     water   low  19       5                 1         6 0.8080
## 620 697     water   low  20       5                 0         0 0.0000
## 621 501     water  high   1      10                 0         0 0.0000
## 622 506     water  high   2      10                 0         0 0.0000
## 623 512     water  high   3      10                 0         0 0.0000
## 624 516     water  high   4      10                 1         4 0.8079
## 625 525     water  high   5      10                 0         0 0.0000
## 626 527     water  high   6      10                 0         0 0.0000
## 627 533     water  high   7      10                 2         9 0.7831
## 628 536     water  high   8      10                 0         0 0.0000
## 629 543     water  high   9      10                 1         7 0.7919
## 630 550     water  high  10      10                 0         0 0.0000
## 631 555     water  high  11      10                 0         0 0.0000
## 632 558     water  high  12      10                 1         5 0.7547
## 633 565     water  high  13      10                 0         0 0.0000
## 634 568     water  high  14      10                 0         0 0.0000
## 635 575     water  high  15      10                 0         0 0.0000
## 636 578     water  high  16      10                 0         0 0.0000
## 637 584     water  high  17      10                 0         0 0.0000
## 638 590     water  high  18      10                 0         0 0.0000
## 639 593     water  high  19      10                 0         0 0.0000
## 640 599     water  high  20      10                 0         0 0.0000
## 641 601     water   low   1      10                 1         7 0.7586
## 642 606     water   low   2      10                 1         5 0.8027
## 643 612     water   low   3      10                 0         0 0.0000
## 644 616     water   low   4      10                 0         0 0.0000
## 645 625     water   low   5      10                 0         0 0.0000
## 646 627     water   low   6      10                 0         0 0.0000
## 647 633     water   low   7      10                 1         6 0.7595
## 648 636     water   low   8      10                 1        10 0.7776
## 649 643     water   low   9      10                 0         0 0.0000
## 650 650     water   low  10      10                 0         0 0.0000
## 651 655     water   low  11      10                 0         0 0.0000
## 652 658     water   low  12      10                 0         0 0.0000
## 653 665     water   low  13      10                 1         7 0.6986
## 654 668     water   low  14      10                 0         0 0.0000
## 655 675     water   low  15      10                 0         0 0.0000
## 656 678     water   low  16      10                 1         7 0.8063
## 657 684     water   low  17      10                 0         0 0.0000
## 658 690     water   low  18      10                 1         7 0.7980
## 659 693     water   low  19      10                 0         0 0.0000
## 660 699     water   low  20      10                 0         0 0.0000
## 661 502     water  high   1      20                 1         7 0.7895
## 662 507     water  high   2      20                 1         7 0.8028
## 663 511     water  high   3      20                 0         0 0.0000
## 664 519     water  high   4      20                 0         0 0.0000
## 665 524     water  high   5      20                 2         8 0.7687
## 666 530     water  high   6      20                 0         0 0.0000
## 667 531     water  high   7      20                 1         5 0.7742
## 668 537     water  high   8      20                 1         7 0.8034
## 669 545     water  high   9      20                 0         0 0.0000
## 670 546     water  high  10      20                 0         0 0.0000
## 671 551     water  high  11      20                 0         0 0.0000
## 672 557     water  high  12      20                 1         2 0.7315
## 673 561     water  high  13      20                 0         0 0.0000
## 674 567     water  high  14      20                 0         0 0.0000
## 675 572     water  high  15      20                 0         0 0.0000
## 676 579     water  high  16      20                 0         0 0.0000
## 677 582     water  high  17      20                 0         0 0.0000
## 678 589     water  high  18      20                 0         0 0.0000
## 679 594     water  high  19      20                 0         0 0.0000
## 680 600     water  high  20      20                 0         0 0.0000
## 681 602     water   low   1      20                 0         0 0.0000
## 682 607     water   low   2      20                 1         6 0.7758
## 683 611     water   low   3      20                 0         0 0.0000
## 684 619     water   low   4      20                 1         7 0.8032
## 685 624     water   low   5      20                 1         5 0.8062
## 686 630     water   low   6      20                 0         0 0.0000
## 687 631     water   low   7      20                 0         0 0.0000
## 688 637     water   low   8      20                 0         0 0.0000
## 689 645     water   low   9      20                 0         0 0.0000
## 690 646     water   low  10      20                 1         5 0.7975
## 691 651     water   low  11      20                 0         0 0.0000
## 692 657     water   low  12      20                 0         0 0.0000
## 693 661     water   low  13      20                 0         0 0.0000
## 694 667     water   low  14      20                 0         0 0.0000
## 695 672     water   low  15      20                 0         0 0.0000
## 696 679     water   low  16      20                 0         0 0.0000
## 697 682     water   low  17      20                 0         0 0.0000
## 698 689     water   low  18      20                 0         0 0.0000
## 699 694     water   low  19      20                 0         0 0.0000
## 700 700     water   low  20      20                 0         0 0.0000
##     brome.begin brome.prop.begin ephedra.begin ephedra.end brome.respawn
## 1             0                0             0           0             0
## 2             0                0             1           1             0
## 3             0                0             1           1             0
## 4             0                0             0           0             0
## 5             0                0             1           0             0
## 6             0                0             0           0             0
## 7             0                0             0           0             0
## 8             0                0             1           1             2
## 9             0                0             0           0             0
## 10            0                0             0           0             0
## 11            0                0             1           1             0
## 12            0                0             0           0             0
## 13            0                0             0           0             0
## 14            0                0             0           0             0
## 15            0                0             1           0             0
## 16            0                0             1           0             0
## 17            0                0             0           0             0
## 18            0                0             1           1             0
## 19            0                0             0           0             0
## 20            0                0             0           0             0
## 21            0                0             0           0             0
## 22            0                0             1           1             0
## 23            0                0             1           1             0
## 24            0                0             1           1             0
## 25            0                0             1           1             0
## 26            0                0             1           1             0
## 27            0                0             1           1             0
## 28            0                0             1           1             0
## 29            0                0             1           1             0
## 30            0                0             1           1             0
## 31            0                0             1           1             0
## 32            0                0             0           0             0
## 33            0                0             1           1             0
## 34            0                0             0           0             0
## 35            0                0             0           0             0
## 36            0                0             1           1             0
## 37            0                0             0           0             0
## 38            0                0             0           0             0
## 39            0                0             0           0             0
## 40            0                0             0           0             0
## 41            1               50             0           0             1
## 42            2              100             0           0             0
## 43            1               50             0           0             1
## 44            2              100             0           0             2
## 45            0                0             0           0             0
## 46            1               50             0           0             1
## 47            2              100             0           0             1
## 48            2              100             0           0             0
## 49            1               50             1           0             0
## 50            1               50             0           0             0
## 51            1               50             0           0             0
## 52            2              100             0           0             2
## 53            2              100             0           0             2
## 54            1               50             1           1             2
## 55            1               50             1           1             0
## 56            2              100             1           1             0
## 57            2              100             0           0             0
## 58            0                0             1           0             0
## 59            0                0             1           0             0
## 60            0                0             0           0             0
## 61            0                0             1           1             0
## 62            1               50             0           0             0
## 63            2              100             1           1             1
## 64            2              100             0           0             1
## 65            2              100             0           0             0
## 66            1               50             0           0             0
## 67            0                0             1           1             1
## 68            2              100             0           0             0
## 69            2              100             0           0             1
## 70            1               50             0           0             1
## 71            2              100             1           1             0
## 72            2              100             0           0             1
## 73            0                0             1           1             0
## 74            2              100             1           1             2
## 75            2              100             0           0             0
## 76            1               50             1           1             0
## 77            1               50             0           0             0
## 78            2              100             0           0             2
## 79            2              100             0           0             0
## 80            0                0             0           0             0
## 81            4               80             1           1             0
## 82            5              100             0           0             0
## 83            5              100             1           1             1
## 84            5              100             0           0             1
## 85            3               60             1           1             0
## 86            4               80             1           1             0
## 87            2               40             1           1             0
## 88            5              100             1           1             1
## 89            4               80             1           1             2
## 90            3               60             1           1             0
## 91            2               40             1           1             1
## 92            1               20             1           1             1
## 93            3               60             1           1             1
## 94            4               80             1           1             0
## 95            2               40             1           1             0
## 96            5              100             0           0             0
## 97            1               20             0           0             2
## 98            5              100             1           1             1
## 99            3               60             1           0             0
## 100           4               80             1           1             1
## 101           5              100             1           1             3
## 102           5              100             1           1             0
## 103           4               80             0           0             2
## 104           3               60             0           0             0
## 105           4               80             1           1             0
## 106           5              100             1           1             1
## 107           2               40             1           1             1
## 108           5              100             1           0             0
## 109           2               40             0           0             0
## 110           4               80             0           0             0
## 111           3               60             0           0             0
## 112           3               60             0           0             2
## 113           0                0             1           1             0
## 114           1               20             0           0             2
## 115           3               60             0           0             3
## 116           5              100             0           0             0
## 117           5              100             1           1             2
## 118           5              100             0           0             0
## 119           4               80             1           1             2
## 120           5              100             0           0             1
## 121           3               30             1           1             0
## 122           5               50             1           1             0
## 123           8               80             0           0             0
## 124           4               40             0           0             0
## 125           8               80             0           0             0
## 126           9               90             0           0             0
## 127          10              100             1           1             2
## 128           5               50             1           1             0
## 129           7               70             0           0             0
## 130           9               90             1           0             1
## 131           3               30             1           1             0
## 132           7               70             1           0             2
## 133           4               40             0           0             0
## 134           3               30             0           0             0
## 135           9               90             1           1             3
## 136           7               70             1           1             0
## 137           6               60             0           0             0
## 138           7               70             1           0             0
## 139           4               40             0           0             1
## 140           0                0             0           0             0
## 141           5               50             0           0             2
## 142           7               70             1           1             0
## 143           7               70             0           0             2
## 144           6               60             0           0             3
## 145           9               90             0           0             3
## 146           6               60             1           1             0
## 147           8               80             0           0             0
## 148           3               30             0           0             1
## 149           8               80             0           0             0
## 150           9               90             1           1             1
## 151           9               90             0           0             0
## 152           5               50             1           1             1
## 153           7               70             1           0             1
## 154          10              100             1           1             0
## 155           8               80             0           0             2
## 156           6               60             1           1             1
## 157           7               70             0           0             0
## 158          10              100             1           1             3
## 159           7               70             0           0             0
## 160           3               30             0           0             2
## 161           8               40             1           1             2
## 162          11               55             0           0             0
## 163          11               55             1           0             3
## 164           9               45             0           0             1
## 165          11               55             0           0             4
## 166           9               45             0           0             3
## 167          13               65             1           0             2
## 168          11               55             1           1             0
## 169           0                0             0           0             0
## 170          13               65             0           0             0
## 171          10               50             0           0             0
## 172           8               40             0           0             0
## 173          10               50             1           1             0
## 174           8               40             0           0             0
## 175           8               40             0           0             0
## 176           9               45             1           0             2
## 177          11               55             0           0             1
## 178          14               70             0           0             2
## 179          10               50             0           0             4
## 180          14               70             1           1             2
## 181           9               45             0           0             2
## 182           9               45             0           0             0
## 183          15               75             0           0             3
## 184          14               70             0           0             3
## 185           9               45             1           1             1
## 186          13               65             1           1             1
## 187          12               60             0           0             2
## 188           1                5             0           0             0
## 189          17               85             0           0             3
## 190          12               60             0           0             3
## 191           8               40             1           0             2
## 192          11               55             1           1             1
## 193          16               80             1           0             0
## 194          10               50             1           1             3
## 195          14               70             0           0             0
## 196          11               55             0           0             2
## 197          10               50             1           1             0
## 198          12               60             0           0             2
## 199           9               45             0           0             4
## 200          13               65             0           0             3
## 201           0                0             1           0             0
## 202           0                0             1           1             0
## 203           0                0             1           1             0
## 204           0                0             1           1             0
## 205           0                0             0           0             0
## 206           0                0             0           0             0
## 207           0                0             1           1             0
## 208           0                0             0           0             0
## 209           0                0             1           1             0
## 210           0                0             1           1             0
## 211           0                0             1           1             0
## 212           0                0             0           0             0
## 213           0                0             1           1             1
## 214           0                0             0           0             0
## 215           0                0             1           0             0
## 216           0                0             1           1             0
## 217           0                0             1           0             0
## 218           0                0             0           0             0
## 219           0                0             0           0             0
## 220           0                0             1           0             0
## 221           2              100             1           1             0
## 222           2              100             0           0             0
## 223           2              100             1           0             0
## 224           2              100             0           0             0
## 225           1               50             0           0             0
## 226           1               50             0           0             1
## 227           2              100             1           1             2
## 228           2              100             1           1             0
## 229           1               50             1           1             0
## 230           1               50             1           1             0
## 231           2              100             1           1             0
## 232           2              100             1           1             0
## 233           1               50             0           0             0
## 234           2              100             0           0             0
## 235           0                0             0           0             0
## 236           2              100             0           0             0
## 237           1               50             1           1             0
## 238           1               50             0           0             0
## 239           1               50             1           0             1
## 240           1               50             0           0             0
## 241           4               80             0           0             0
## 242           3               60             0           0             3
## 243           5              100             0           0             0
## 244           4               80             0           0             1
## 245           3               60             0           0             1
## 246           2               40             1           1             1
## 247           4               80             1           1             1
## 248           3               60             1           0             3
## 249           3               60             1           1             2
## 250           2               40             1           0             0
## 251           5              100             0           0             2
## 252           3               60             0           0             0
## 253           5              100             1           1             0
## 254           3               60             1           1             1
## 255           4               80             1           0             0
## 256           3               60             1           1             0
## 257           4               80             1           0             0
## 258           3               60             0           0             0
## 259           5              100             1           1             2
## 260           5              100             0           0             1
## 261           9               90             0           0             1
## 262           5               50             1           1             1
## 263           5               50             0           0             2
## 264           6               60             0           0             0
## 265           5               50             1           0             0
## 266          10              100             1           1             2
## 267           7               70             1           1             0
## 268           4               40             0           0             0
## 269           8               80             1           0             0
## 270           5               50             0           0             0
## 271           7               70             1           1             0
## 272           3               30             0           0             0
## 273           7               70             1           0             0
## 274           8               80             0           0             0
## 275           5               50             1           1             0
## 276           7               70             0           0             0
## 277           8               80             1           1             1
## 278          10              100             0           0             1
## 279           4               40             1           0             0
## 280           7               70             0           0             1
## 281          10               50             1           1             2
## 282          10               50             1           1             2
## 283          14               70             0           0             3
## 284          10               50             0           0             1
## 285          12               60             0           0             2
## 286          10               50             0           0             2
## 287          12               60             1           1             0
## 288          10               50             1           1             0
## 289           7               35             0           0             1
## 290          10               50             0           0             5
## 291          14               70             1           0             2
## 292           6               30             1           1             0
## 293           9               45             0           0             0
## 294          14               70             0           0             0
## 295           9               45             1           0             0
## 296          12               60             0           0             2
## 297          15               75             1           0             0
## 298          12               60             1           0             1
## 299           8               40             1           0             2
## 300          13               65             1           1             0
## 301           0                0             0           0             0
## 302           0                0             1           1             0
## 303           0                0             0           0             0
## 304           0                0             1           1             0
## 305           0                0             0           0             0
## 306           0                0             1           1             0
## 307           0                0             0           0             0
## 308           0                0             0           0             0
## 309           0                0             1           1             0
## 310           0                0             0           0             0
## 311           0                0             0           0             0
## 312           0                0             0           0             0
## 313           0                0             1           1             0
## 314           0                0             0           0             0
## 315           0                0             0           0             0
## 316           0                0             0           0             0
## 317           0                0             1           1             0
## 318           0                0             0           0             0
## 319           0                0             1           0             0
## 320           0                0             1           1             0
## 321           0                0             0           0             0
## 322           0                0             0           0             0
## 323           0                0             0           0             0
## 324           0                0             0           0             0
## 325           0                0             0           0             0
## 326           0                0             1           0             0
## 327           0                0             1           1             0
## 328           0                0             1           0             0
## 329           0                0             0           0             0
## 330           0                0             1           0             0
## 331           0                0             0           0             0
## 332           0                0             0           0             0
## 333           0                0             0           0             0
## 334           0                0             0           0             0
## 335           0                0             0           0             0
## 336           0                0             1           0             0
## 337           0                0             0           0             0
## 338           0                0             1           1             0
## 339           0                0             1           0             0
## 340           0                0             0           0             0
## 341           2              100             0           0             0
## 342           1               50             0           0             0
## 343           2              100             0           0             0
## 344           2              100             1           1             0
## 345           2              100             1           1             0
## 346           2              100             0           0             0
## 347           2              100             0           0             0
## 348           2              100             1           1             0
## 349           1               50             0           0             0
## 350           2              100             0           0             2
## 351           2              100             0           0             0
## 352           2              100             0           0             0
## 353           2              100             1           1             0
## 354           2              100             0           0             0
## 355           2              100             1           1             1
## 356           2              100             0           0             0
## 357           2              100             0           0             0
## 358           2              100             1           1             0
## 359           2              100             0           0             0
## 360           2              100             0           0             0
## 361           1               50             1           0             0
## 362           1               50             0           0             0
## 363           2              100             0           0             0
## 364           0                0             0           0             0
## 365           1               50             1           0             0
## 366           1               50             1           1             0
## 367           0                0             0           0             0
## 368           2              100             0           0             0
## 369           1               50             1           1             0
## 370           2              100             1           0             0
## 371           1               50             0           0             0
## 372           2              100             1           1             0
## 373           2              100             0           0             0
## 374           2              100             0           0             0
## 375           1               50             1           1             0
## 376           2              100             0           0             0
## 377           1               50             1           1             0
## 378           0                0             1           1             0
## 379           0                0             0           0             0
## 380           0                0             0           0             0
## 381           4               80             1           1             0
## 382           5              100             1           1             0
## 383           4               80             1           1             0
## 384           3               60             1           1             0
## 385           4               80             1           1             0
## 386           4               80             0           0             0
## 387           4               80             0           0             0
## 388           4               80             0           0             0
## 389           3               60             0           0             0
## 390           5              100             1           1             0
## 391           4               80             0           0             0
## 392           5              100             0           0             0
## 393           2               40             1           1             0
## 394           4               80             0           0             2
## 395           4               80             0           0             0
## 396           4               80             0           0             0
## 397           5              100             0           0             0
## 398           3               60             1           1             0
## 399           3               60             1           1             0
## 400           4               80             0           0             0
## 401           3               60             0           0             0
## 402           3               60             1           1             0
## 403           1               20             1           1             0
## 404           4               80             1           1             0
## 405           3               60             0           0             0
## 406           5              100             1           0             0
## 407           5              100             0           0             3
## 408           2               40             1           1             0
## 409           3               60             0           0             0
## 410           3               60             0           0             2
## 411           4               80             0           0             0
## 412           4               80             0           0             0
## 413           4               80             1           1             0
## 414           3               60             0           0             0
## 415           0                0             1           0             0
## 416           1               20             0           0             0
## 417           4               80             0           0             0
## 418           3               60             0           0             0
## 419           4               80             0           0             0
## 420           0                0             0           0             0
## 421          10              100             1           1             1
## 422           7               70             0           0             0
## 423           7               70             0           0             0
## 424           7               70             0           0             0
## 425           3               30             0           0             0
## 426           7               70             0           0             0
## 427           8               80             1           1             0
## 428           7               70             1           1             0
## 429           8               80             0           0             0
## 430           7               70             0           0             0
## 431           5               50             0           0             0
## 432           5               50             1           1             0
## 433           7               70             0           0             0
## 434           6               60             0           0             0
## 435           5               50             0           0             0
## 436           8               80             1           1             0
## 437           7               70             0           0             0
## 438           4               40             0           0             0
## 439           5               50             1           1             0
## 440           9               90             1           1             0
## 441           5               50             1           1             0
## 442          10              100             0           0             1
## 443           6               60             0           0             0
## 444           8               80             1           1             0
## 445           6               60             0           0             0
## 446           6               60             0           0             0
## 447           6               60             0           0             0
## 448           7               70             0           0             1
## 449           7               70             0           0             1
## 450           6               60             0           0             0
## 451           5               50             0           0             1
## 452           5               50             0           0             0
## 453           7               70             0           0             0
## 454           7               70             1           1             0
## 455           7               70             0           0             0
## 456           4               40             0           0             0
## 457           4               40             0           0             0
## 458           8               80             0           0             0
## 459           5               50             0           0             0
## 460           5               50             0           0             0
## 461           8               40             0           0             0
## 462          14               70             1           1             0
## 463          11               55             0           0             0
## 464          12               60             0           0             0
## 465          14               70             1           1             0
## 466           9               45             0           0             0
## 467          11               55             1           1             0
## 468           8               40             0           0             0
## 469          12               60             0           0             0
## 470          14               70             1           1             2
## 471          12               60             1           1             3
## 472          13               65             0           0             3
## 473           8               40             0           0             0
## 474          12               60             1           1             0
## 475          10               50             0           0             0
## 476          16               80             1           1             0
## 477          13               65             0           0             1
## 478          11               55             0           0             0
## 479          12               60             1           1             0
## 480           9               45             0           0             0
## 481           8               40             0           0             0
## 482          12               60             0           0             0
## 483           9               45             0           0             0
## 484          12               60             0           0             1
## 485           8               40             0           0             0
## 486          13               65             0           0             1
## 487           9               45             0           0             0
## 488           7               35             0           0             1
## 489           7               35             0           0             0
## 490           3               15             0           0             0
## 491           6               30             1           1             0
## 492           8               40             0           0             0
## 493          10               50             1           1             1
## 494          10               50             0           0             0
## 495          10               50             0           0             1
## 496           1                5             0           0             0
## 497           9               45             0           0             0
## 498          10               50             0           0             0
## 499          10               50             0           0             0
## 500           7               35             0           0             1
## 501           0                0             1           1             0
## 502           0                0             1           1             0
## 503           0                0             1           1             0
## 504           0                0             1           1             0
## 505           0                0             1           1             0
## 506           0                0             0           0             0
## 507           0                0             0           0             0
## 508           0                0             0           0             0
## 509           0                0             1           0             0
## 510           0                0             1           1             0
## 511           0                0             1           0             0
## 512           0                0             0           0             0
## 513           0                0             1           1             0
## 514           0                0             1           0             0
## 515           0                0             1           1             0
## 516           0                0             1           0             0
## 517           0                0             0           0             0
## 518           0                0             1           0             0
## 519           0                0             1           1             0
## 520           0                0             1           0             0
## 521           0                0             1           1             0
## 522           0                0             1           1             0
## 523           0                0             1           1             0
## 524           0                0             1           1             0
## 525           0                0             0           0             0
## 526           0                0             1           1             0
## 527           0                0             1           1             0
## 528           0                0             0           0             0
## 529           0                0             1           0             0
## 530           0                0             1           1             0
## 531           0                0             1           0             0
## 532           0                0             0           0             0
## 533           0                0             0           0             0
## 534           0                0             1           1             0
## 535           0                0             0           0             0
## 536           0                0             0           0             0
## 537           0                0             0           0             0
## 538           0                0             1           1             0
## 539           0                0             0           0             0
## 540           0                0             0           0             0
## 541           1               50             1           0             0
## 542           2              100             0           0             0
## 543           2              100             1           0             2
## 544           1               50             0           0             0
## 545           1               50             1           1             0
## 546           0                0             0           0             2
## 547           2              100             1           1             0
## 548           1               50             1           1             0
## 549           0                0             1           1             0
## 550           1               50             1           0             0
## 551           1               50             1           1             0
## 552           2              100             1           0             0
## 553           2              100             0           0             0
## 554           2              100             1           0             0
## 555           0                0             1           1             0
## 556           2              100             1           1             0
## 557           1               50             1           1             0
## 558           2              100             1           1             0
## 559           2              100             1           1             0
## 560           2              100             1           0             0
## 561           2              100             1           1             0
## 562           2              100             1           1             0
## 563           2              100             0           0             0
## 564           2              100             0           0             0
## 565           1               50             0           0             0
## 566           2              100             1           1             0
## 567           2              100             0           0             0
## 568           2              100             0           0             0
## 569           1               50             1           1             0
## 570           0                0             0           0             0
## 571           0                0             1           0             0
## 572           1               50             0           0             0
## 573           1               50             1           1             0
## 574           2              100             1           1             0
## 575           2              100             1           0             0
## 576           2              100             0           0             0
## 577           2              100             1           1             0
## 578           2              100             0           0             0
## 579           1               50             0           0             0
## 580           2              100             1           1             0
## 581           4               80             1           1             0
## 582           0                0             1           1             0
## 583           5              100             1           1             0
## 584           3               60             1           1             0
## 585           3               60             1           1             0
## 586           4               80             0           0             0
## 587           0                0             0           0             0
## 588           5              100             0           0             3
## 589           5              100             0           0             0
## 590           3               60             0           0             0
## 591           5              100             1           1             0
## 592           2               40             0           0             0
## 593           5              100             0           0             0
## 594           3               60             0           0             0
## 595           2               40             0           0             0
## 596           1               20             1           1             0
## 597           5              100             1           1             0
## 598           2               40             1           0             2
## 599           5              100             0           0             1
## 600           5              100             0           0             0
## 601           4               80             0           0             0
## 602           5              100             0           0             0
## 603           5              100             0           0             0
## 604           5              100             0           0             1
## 605           4               80             1           1             0
## 606           5              100             1           1             0
## 607           3               60             0           0             0
## 608           5              100             0           0             0
## 609           3               60             0           0             0
## 610           4               80             0           0             0
## 611           5              100             0           0             2
## 612           3               60             0           0             0
## 613           4               80             0           0             1
## 614           3               60             0           0             0
## 615           4               80             1           0             0
## 616           2               40             0           0             0
## 617           5              100             1           1             0
## 618           3               60             1           1             0
## 619           3               60             1           1             0
## 620           0                0             1           1             0
## 621           2               20             0           0             0
## 622           6               60             1           0             0
## 623           6               60             0           0             0
## 624           6               60             1           1             0
## 625           6               60             0           0             0
## 626           5               50             0           0             0
## 627           7               70             1           1             0
## 628           5               50             0           0             0
## 629           8               80             1           1             0
## 630           6               60             0           0             0
## 631           8               80             1           0             0
## 632           5               50             1           1             0
## 633           5               50             0           0             0
## 634           6               60             0           0             0
## 635           9               90             0           0             3
## 636           6               60             0           0             0
## 637           1               10             1           1             0
## 638           4               40             1           0             2
## 639           5               50             1           0             0
## 640           5               50             0           0             0
## 641           7               70             1           1             0
## 642           4               40             1           1             0
## 643           7               70             1           0             0
## 644           7               70             0           0             0
## 645           8               80             1           1             0
## 646          10              100             0           0             0
## 647           8               80             1           1             0
## 648           5               50             1           1             0
## 649           6               60             0           0             0
## 650          10              100             0           0             0
## 651           7               70             0           0             0
## 652           4               40             0           0             0
## 653           7               70             1           1             0
## 654           6               60             1           1             0
## 655           3               30             0           0             0
## 656           7               70             1           1             0
## 657           7               70             0           0             0
## 658           4               40             1           1             1
## 659           7               70             0           0             1
## 660           7               70             0           0             0
## 661          12               60             1           0             0
## 662          13               65             1           1             0
## 663          13               65             1           1             0
## 664           8               40             1           1             0
## 665          15               75             1           1             0
## 666           5               25             1           0             0
## 667          13               65             1           0             1
## 668           7               35             1           1             0
## 669           8               40             0           0             0
## 670           7               35             1           0             0
## 671           6               30             0           0             0
## 672           1                5             1           1             0
## 673           8               40             1           0             0
## 674          13               65             0           0             0
## 675          11               55             1           0             0
## 676          11               55             0           0             0
## 677          10               50             1           1             0
## 678          13               65             0           0             4
## 679          11               55             0           0             0
## 680          12               60             0           0             0
## 681           7               35             0           0             0
## 682          10               50             1           1             0
## 683          15               75             0           0             0
## 684          14               70             1           1             0
## 685          10               50             1           1             1
## 686          15               75             1           0             1
## 687           8               40             0           0             0
## 688          13               65             0           0             0
## 689          13               65             1           0             0
## 690          12               60             1           0             0
## 691          10               50             0           0             0
## 692           7               35             0           0             0
## 693          12               60             1           1             2
## 694           8               40             0           0             0
## 695           9               45             0           0             0
## 696          13               65             0           0             1
## 697          11               55             1           0             0
## 698          13               65             1           1             2
## 699          14               70             0           0             0
## 700          10               50             1           1             0
##     Ephedra.above Ephedra.below Ephedra.biomass brome.biomass
## 1         0.00550       0.00120         0.00670       0.00000
## 2         0.02650       0.00350         0.03000       0.09180
## 3         0.01360       0.00510         0.01870       0.25030
## 4         0.00000       0.00000         0.00000       0.00000
## 5         0.02130       0.00530         0.02660       0.12160
## 6         0.00000       0.00000         0.00000       0.29580
## 7         0.00000       0.00000         0.00000       0.26070
## 8         0.01650       0.00820         0.02470       0.18090
## 9         0.00000       0.00000         0.00000       0.20730
## 10        0.00000       0.00000         0.00000       0.28800
## 11        0.00000       0.00000         0.00000       0.00000
## 12        0.03620       0.00600         0.04220       0.37190
## 13        0.00000       0.00000         0.00000       0.05550
## 14        0.00000       0.00000         0.00000       0.00000
## 15        0.00000       0.00000         0.00000       0.19680
## 16        0.00000       0.00000         0.00000       0.00000
## 17        0.02370       0.00460         0.02830       0.09410
## 18        0.00000       0.00000         0.00000       0.64990
## 19        0.00000       0.00000         0.00000       0.37030
## 20        0.00000       0.00000         0.00000       0.38870
## 21        0.00000       0.00000         0.00000       0.10390
## 22        0.00000       0.00000         0.00000       0.12920
## 23        0.01850       0.00000         0.01850       0.08190
## 24        0.01300       0.00080         0.01380       0.00000
## 25        0.00000       0.00000         0.00000       0.20410
## 26        0.00000       0.00000         0.00000       0.12900
## 27        0.00000       0.00000         0.00000       0.11510
## 28        0.01280       0.00530         0.01810       0.07580
## 29        0.00000       0.00000         0.00000       0.10430
## 30        0.02450       0.00260         0.02710       0.30570
## 31        0.03030       0.00830         0.03860       0.00000
## 32        0.04260       0.00930         0.05190       0.22820
## 33        0.00000       0.00000         0.00000       0.16880
## 34        0.00930       0.00320         0.01250       0.00000
## 35        0.00000       0.00000         0.00000       0.63810
## 36        0.01810       0.00310         0.02120       0.00000
## 37        0.01650       0.00270         0.01920       0.05580
## 38        0.01550       0.00080         0.01630       0.25570
## 39        0.00000       0.00000         0.00000       0.31530
## 40        0.00000       0.00000         0.00000       0.43460
## 41        0.01070       0.00080         0.01150       0.28450
## 42        0.00000       0.00000         0.00000       0.18080
## 43        0.00000       0.00000         0.00000       0.00000
## 44        0.00000       0.00000         0.00000       0.51090
## 45        0.00000       0.00000         0.00000       0.49770
## 46        0.00000       0.00000         0.00000       0.37290
## 47        0.00000       0.00000         0.00000       0.00000
## 48        0.02490       0.00640         0.03130       0.48910
## 49        0.00000       0.00000         0.00000       0.26980
## 50        0.03850       0.00740         0.04590       0.09020
## 51        0.00000       0.00000         0.00000       0.06860
## 52        0.00000       0.00000         0.00000       0.00000
## 53        0.00000       0.00000         0.00000       0.04360
## 54        0.00000       0.00000         0.00000       0.04860
## 55        0.00000       0.00000         0.00000       0.34990
## 56        0.00000       0.00000         0.00000       0.08430
## 57        0.00000       0.00000         0.00000       0.05540
## 58        0.01160       0.00350         0.01510       0.42220
## 59        0.05680       0.00790         0.06470       0.08060
## 60        0.00000       0.00000         0.00000       0.14490
## 61        0.04010       0.01200         0.05210       0.22380
## 62        0.02610       0.01290         0.03900       0.07840
## 63        0.00000       0.00000         0.00000       0.00000
## 64        0.01410       0.00560         0.01970       0.10500
## 65        0.01290       0.00410         0.01700       0.31520
## 66        0.03570       0.00750         0.04320       0.13820
## 67        0.00000       0.00000         0.00000       0.00000
## 68        0.00000       0.00000         0.00000       0.08840
## 69        0.00000       0.00000         0.00000       0.20450
## 70        0.00900       0.00360         0.01260       0.05780
## 71        0.00000       0.00000         0.00000       0.07440
## 72        0.00900       0.00250         0.01150       0.32930
## 73        0.00000       0.00000         0.00000       0.03150
## 74        0.00990       0.00500         0.01490       0.00000
## 75        0.00000       0.00000         0.00000       0.29070
## 76        0.00000       0.00000         0.00000       0.06400
## 77        0.00000       0.00000         0.00000       0.36160
## 78        0.00000       0.00000         0.00000       0.08940
## 79        0.02130       0.00510         0.02640       0.06920
## 80        0.00000       0.00000         0.00000       0.14790
## 81        0.00000       0.00000         0.00000       0.97380
## 82        0.00000       0.00000         0.00000       0.07510
## 83        0.00000       0.00000         0.00000       0.00000
## 84        0.02700       0.00230         0.02930       0.37440
## 85        0.00000       0.00000         0.00000       0.00000
## 86        0.03300       0.00550         0.03850       0.05470
## 87        0.00000       0.00000         0.00000       0.28710
## 88        0.00000       0.00000         0.00000       0.13070
## 89        0.00000       0.00000         0.00000       0.00000
## 90        0.00000       0.00000         0.00000       0.00000
## 91        0.04280       0.00770         0.05050       0.14800
## 92        0.00000       0.00000         0.00000       0.00000
## 93        0.00000       0.00000         0.00000       0.06470
## 94        0.00000       0.00000         0.00000       0.36080
## 95        0.00000       0.00000         0.00000       0.00000
## 96        0.00560       0.00270         0.00830       0.48060
## 97        0.00000       0.00000         0.00000       0.18510
## 98        0.00000       0.00000         0.00000       0.03410
## 99        0.00850       0.00180         0.01030       0.24800
## 100       0.01120       0.00560         0.01680       0.12510
## 101       0.00000       0.00000         0.00000       0.22280
## 102       0.00000       0.00000         0.00000       0.27740
## 103       0.00000       0.00000         0.00000       0.00000
## 104       0.00000       0.00000         0.00000       0.16200
## 105       0.01540       0.00220         0.01760       0.00000
## 106       0.00000       0.00000         0.00000       0.13150
## 107       0.01820       0.00520         0.02340       0.13440
## 108       0.04970       0.01330         0.06300       0.19300
## 109       0.01480       0.00230         0.01710       0.00000
## 110       0.01200       0.00220         0.01420       0.04960
## 111       0.00000       0.00000         0.00000       0.15610
## 112       0.00000       0.00000         0.00000       0.49510
## 113       0.00820       0.00220         0.01040       0.07590
## 114       0.00000       0.00000         0.00000       0.55020
## 115       0.00000       0.00000         0.00000       0.00000
## 116       0.00000       0.00000         0.00000       0.74260
## 117       0.00000       0.00000         0.00000       0.23590
## 118       0.00000       0.00000         0.00000       0.06120
## 119       0.00150       0.00160         0.00310       0.09730
## 120       0.00000       0.00000         0.00000       0.07230
## 121       0.02430       0.00680         0.03110       0.34550
## 122       0.00000       0.00000         0.00000       0.08370
## 123       0.00000       0.00000         0.00000       0.52100
## 124       0.02340       0.00480         0.02820       0.17470
## 125       0.00000       0.00000         0.00000       0.00000
## 126       0.00000       0.00000         0.00000       0.36200
## 127       0.00000       0.00000         0.00000       0.00000
## 128       0.00000       0.00000         0.00000       0.11290
## 129       0.02560       0.00530         0.03090       0.11980
## 130       0.01280       0.00260         0.01540       0.29680
## 131       0.00000       0.00000         0.00000       0.03630
## 132       0.01900       0.00490         0.02390       0.12700
## 133       0.02270       0.00920         0.03190       0.03560
## 134       0.04150       0.00480         0.04630       0.31550
## 135       0.00000       0.00000         0.00000       0.00000
## 136       0.00000       0.00000         0.00000       0.00000
## 137       0.00000       0.00000         0.00000       0.10080
## 138       0.00000       0.00000         0.00000       0.00000
## 139       0.00000       0.00000         0.00000       0.75300
## 140       0.00000       0.00000         0.00000       0.11050
## 141       0.00000       0.00000         0.00000       0.11840
## 142       0.02530       0.00610         0.03140       0.07840
## 143       0.00000       0.00000         0.00000       0.04040
## 144       0.00000       0.00000         0.00000       0.21310
## 145       0.00820       0.00490         0.01310       0.00000
## 146       0.00300       0.00060         0.00360       0.20720
## 147       0.00000       0.00000         0.00000       0.11400
## 148       0.00000       0.00000         0.00000       0.06970
## 149       0.00000       0.00000         0.00000       0.08230
## 150       0.00000       0.00000         0.00000       0.26620
## 151       0.00000       0.00000         0.00000       0.05840
## 152       0.00000       0.00000         0.00000       0.04610
## 153       0.00000       0.00000         0.00000       0.09150
## 154       0.00000       0.00000         0.00000       0.50060
## 155       0.00000       0.00000         0.00000       0.00000
## 156       0.00000       0.00000         0.00000       0.00000
## 157       0.00000       0.00000         0.00000       0.11740
## 158       0.01170       0.00270         0.01440       0.00000
## 159       0.00000       0.00000         0.00000       0.40220
## 160       0.01240       0.00200         0.01440       0.05850
## 161       0.00000       0.00000         0.00000       0.75040
## 162       0.06560       0.01170         0.07730       0.00000
## 163       0.00000       0.00000         0.00000       0.00000
## 164       0.03600       0.00780         0.04380       0.13120
## 165       0.00000       0.00000         0.00000       0.33800
## 166       0.01150       0.00230         0.01380       0.12240
## 167       0.00000       0.00000         0.00000       0.22720
## 168       0.00000       0.00000         0.00000       0.18500
## 169       0.03550       0.00100         0.03650       0.14660
## 170       0.01950       0.00380         0.02330       0.29450
## 171       0.00810       0.00500         0.01310       0.40970
## 172       0.01220       0.00430         0.01650       0.00000
## 173       0.00000       0.00000         0.00000       0.35110
## 174       0.00000       0.00000         0.00000       0.59720
## 175       0.00000       0.00000         0.00000       0.00000
## 176       0.00000       0.00000         0.00000       0.00000
## 177       0.00000       0.00000         0.00000       0.20830
## 178       0.00000       0.00000         0.00000       0.00000
## 179       0.00000       0.00000         0.00000       0.00000
## 180       0.00000       0.00000         0.00000       0.14880
## 181       0.00770       0.00000         0.00770       0.16940
## 182       0.01310       0.00270         0.01580       0.00000
## 183       0.00000       0.00000         0.00000       0.00000
## 184       0.00000       0.00000         0.00000       0.08240
## 185       0.00000       0.00000         0.00000       0.22150
## 186       0.00720       0.00080         0.00800       0.08980
## 187       0.00000       0.00000         0.00000       0.05690
## 188       0.00000       0.00000         0.00000       0.05500
## 189       0.00000       0.00000         0.00000       0.63270
## 190       0.00000       0.00000         0.00000       0.30780
## 191       0.03100       0.00870         0.03970       0.23010
## 192       0.01410       0.00130         0.01540       0.00000
## 193       0.01750       0.00130         0.01880       0.46290
## 194       0.01240       0.00260         0.01500       0.21450
## 195       0.00000       0.00000         0.00000       0.00000
## 196       0.00000       0.00000         0.00000       0.00000
## 197       0.00000       0.00000         0.00000       0.52390
## 198       0.00000       0.00000         0.00000       0.00000
## 199       0.00000       0.00000         0.00000       0.05120
## 200       0.01320       0.00140         0.01460       0.00000
## 201       0.00000       0.00000         0.00000       0.26740
## 202       0.01110       0.00190         0.01300       0.15880
## 203       0.03880       0.00270         0.04150       0.10360
## 204       0.00000       0.00000         0.00000       0.00000
## 205       0.00000       0.00000         0.00000       0.00000
## 206       0.03070       0.00230         0.03300       0.12130
## 207       0.00000       0.00000         0.00000       0.19280
## 208       0.00000       0.00000         0.00000       0.17980
## 209       0.69900       0.01180         0.07180       0.79270
## 210       0.03840       0.00690         0.04530       0.58480
## 211       0.00000       0.00000         0.00000       0.00000
## 212       0.00000       0.00000         0.00000       0.54080
## 213       0.00000       0.00000         0.00000       0.28460
## 214       0.04600       0.00760         0.05360       0.00000
## 215       0.00000       0.00000         0.00000       0.64600
## 216       0.02830       0.00550         0.03380       0.00000
## 217       0.00000       0.00000         0.00000       0.18060
## 218       0.00000       0.00000         0.00000       0.51500
## 219       0.00000       0.00000         0.00000       0.47590
## 220       0.00000       0.00000         0.00000       0.62080
## 221       0.01240       0.00530         0.01770       0.33220
## 222       0.08860       0.01020         0.09880       0.24710
## 223       0.00000       0.00000         0.00000       0.23960
## 224       0.00350       0.00000         0.00350       0.47170
## 225       0.04330       0.00440         0.04770       0.26980
## 226       0.03900       0.00650         0.04550       0.44520
## 227       0.00000       0.00000         0.00000       0.07190
## 228       0.00000       0.00000         0.00000       0.25540
## 229       0.00000       0.00000         0.00000       0.56060
## 230       0.05980       0.00940         0.06920       0.48800
## 231       0.03560       0.00410         0.03970       0.01880
## 232       0.05700       0.01200         0.06900       0.70010
## 233       0.00000       0.00000         0.00000       0.09610
## 234       0.00000       0.00000         0.00000       0.08420
## 235       0.00000       0.00000         0.00000       0.19430
## 236       0.05690       0.01030         0.06720       0.16430
## 237       0.00000       0.00000         0.00000       0.10320
## 238       0.03300       0.00790         0.04090       0.25460
## 239       0.00000       0.00000         0.00000       0.15270
## 240       0.02760       0.00770         0.03530       0.19710
## 241       0.00000       0.00000         0.00000       0.40440
## 242       0.06080       0.01290         0.07370       0.34680
## 243       0.00000       0.00000         0.00000       0.00000
## 244       0.01850       0.00410         0.02260       0.12600
## 245       0.03220       0.00850         0.04070       0.00000
## 246       0.00000       0.00000         0.00000       0.72250
## 247       0.00000       0.00000         0.00000       0.58350
## 248       0.01660       0.00140         0.01800       0.22040
## 249       0.03030       0.00380         0.03410       0.00000
## 250       0.03990       0.00470         0.04460       0.00000
## 251       0.00000       0.00000         0.00000       0.21270
## 252       0.00000       0.00000         0.00000       0.14510
## 253       0.00000       0.00000         0.00000       0.07130
## 254       0.00000       0.00000         0.00000       0.55520
## 255       0.00000       0.00000         0.00000       0.00000
## 256       0.00000       0.00000         0.00000       0.58590
## 257       0.00000       0.00000         0.00000       0.20620
## 258       0.00000       0.00000         0.00000       0.21800
## 259       0.00000       0.00000         0.00000       0.13050
## 260       0.00000       0.00000         0.00000       0.21200
## 261       0.00000       0.00000         0.00000       0.39010
## 262       0.01580       0.00400         0.01980       0.18070
## 263       0.00000       0.00000         0.00000       0.08230
## 264       0.04550       0.00490         0.05040       0.31070
## 265       0.06410       0.01890         0.08300       0.00000
## 266       0.00000       0.00000         0.00000       0.61500
## 267       0.00000       0.00000         0.00000       0.00000
## 268       0.00000       0.00000         0.00000       0.07120
## 269       0.03720       0.00550         0.04270       0.32410
## 270       0.06100       0.01660         0.07760       0.21200
## 271       0.00000       0.00000         0.00000       0.09400
## 272       0.00000       0.00000         0.00000       0.16790
## 273       0.06320       0.00790         0.07110       0.15580
## 274       0.00000       0.00000         0.00000       0.33890
## 275       0.00000       0.00000         0.00000       0.00000
## 276       0.00000       0.00000         0.00000       0.00000
## 277       0.00920       0.00210         0.01130       0.13280
## 278       0.00000       0.00000         0.00000       0.00000
## 279       0.00000       0.00000         0.00000       0.23810
## 280       0.02680       0.00730         0.03410       0.07970
## 281       0.00000       0.00000         0.00000       0.40410
## 282       0.00000       0.00000         0.00000       0.00000
## 283       0.00000       0.00000         0.00000       0.00000
## 284       0.00000       0.00000         0.00000       0.53350
## 285       0.00000       0.00000         0.00000       0.64980
## 286       0.03330       0.00000         0.03330       0.43210
## 287       0.00000       0.00610         0.00610       0.11750
## 288       0.00000       0.00000         0.00000       0.18620
## 289       0.04940       0.01460         0.06400       0.63010
## 290       0.00000       0.00000         0.00000       0.42480
## 291       0.00000       0.00000         0.00000       0.19640
## 292       0.00000       0.00000         0.00000       0.00000
## 293       0.00000       0.00000         0.00000       0.40760
## 294       0.01220       0.00270         0.01490       0.24880
## 295       0.00000       0.00000         0.00000       0.00000
## 296       0.00000       0.00000         0.00000       0.00000
## 297       0.00000       0.00000         0.00000       0.35760
## 298       0.00000       0.00000         0.00000       0.00000
## 299       0.00000       0.00000         0.00000       0.17540
## 300       0.00000       0.00000         0.00000       0.26390
## 301       0.03520       0.00490         0.04010       0.13690
## 302       0.00000       0.00000         0.00000       0.11630
## 303       0.00000       0.00000         0.00000       0.13920
## 304       0.00000       0.00000         0.00000       0.00000
## 305       0.00000       0.00000         0.00000       0.10280
## 306       0.03260       0.00080         0.03340       0.03700
## 307       0.00000       0.00000         0.00000       0.13280
## 308       0.00000       0.00000         0.00000       0.09610
## 309       0.00000       0.00000         0.00000       0.28290
## 310       0.00000       0.00000         0.00000       0.00000
## 311       0.00000       0.00000         0.00000       0.00000
## 312       0.00000       0.00000         0.00000       0.22240
## 313       0.02280       0.00570         0.02850       0.23670
## 314       0.00000       0.00000         0.00000       0.00000
## 315       0.00000       0.00000         0.00000       0.25500
## 316       0.02230       0.00270         0.02500       0.00000
## 317       0.00000       0.00000         0.00000       0.07660
## 318       0.00000       0.00000         0.00000       0.23500
## 319       0.00000       0.00000         0.00000       0.13050
## 320       0.00000       0.00000         0.00000       0.00000
## 321       0.00000       0.00000         0.00000       0.16580
## 322       0.00000       0.00000         0.00000       0.20550
## 323       0.00000       0.00000         0.00000       0.31370
## 324       0.05320       0.00620         0.05940       0.00000
## 325       0.01200       0.00130         0.01330       0.22010
## 326       0.00000       0.00000         0.00000       0.30350
## 327       0.00000       0.00000         0.00000       0.36030
## 328       0.00000       0.00000         0.00000       0.14580
## 329       0.00000       0.00000         0.00000       0.20870
## 330       0.00000       0.00000         0.00000       0.23930
## 331       0.00000       0.00000         0.00000       0.00000
## 332       0.00000       0.00000         0.00000       0.32050
## 333       0.00000       0.00000         0.00000       0.18900
## 334       0.00000       0.00000         0.00000       0.00000
## 335       0.00000       0.00000         0.00000       0.14870
## 336       0.00000       0.00000         0.00000       0.00000
## 337       0.00000       0.00000         0.00000       0.09670
## 338       0.00000       0.00000         0.00000       0.56820
## 339       0.00000       0.00000         0.00000       0.22460
## 340       0.00000       0.00000         0.00000       0.86120
## 341       0.00000       0.00000         0.00000       0.24180
## 342       0.02870       0.00230         0.03100       0.15810
## 343       0.00000       0.00000         0.00000       0.00000
## 344       0.00000       0.00000         0.00000       0.28430
## 345       0.00000       0.00000         0.00000       0.35110
## 346       0.00000       0.00000         0.00000       1.11320
## 347       0.00000       0.00000         0.00000       0.01380
## 348       0.00000       0.00000         0.00000       0.09180
## 349       0.01830       0.00280         0.02110       0.09270
## 350       0.00550       0.00040         0.00590       0.03290
## 351       0.00000       0.00000         0.00000       0.03820
## 352       0.00000       0.00000         0.00000       0.13830
## 353       0.00000       0.00000         0.00000       0.08510
## 354       0.00000       0.00000         0.00000       0.04030
## 355       0.00000       0.00000         0.00000       0.18130
## 356       0.00000       0.00000         0.00000       0.04520
## 357       0.00000       0.00000         0.00000       0.08300
## 358       0.00000       0.00000         0.00000       0.04690
## 359       0.00000       0.00000         0.00000       0.04090
## 360       0.00000       0.00000         0.00000       0.08500
## 361       0.00000       0.00000         0.00000       0.28350
## 362       0.07030       0.00830         0.07860       0.27440
## 363       0.00000       0.00000         0.00000       0.00000
## 364       0.00000       0.00000         0.00000       0.72950
## 365       0.00000       0.00000         0.00000       0.73950
## 366       0.00000       0.00000         0.00000       0.52640
## 367       0.00000       0.00000         0.00000       0.00000
## 368       0.00000       0.00000         0.00000       0.58350
## 369       0.00000       0.00000         0.00000       0.37550
## 370       0.03490       0.00670         0.04160       0.14020
## 371       0.00000       0.00000         0.00000       0.21470
## 372       0.00000       0.00000         0.00000       0.59210
## 373       0.00000       0.00000         0.00000       0.00000
## 374       0.00000       0.00000         0.00000       0.02350
## 375       0.00000       0.00000         0.00000       0.48560
## 376       0.00000       0.00000         0.00000       0.11270
## 377       0.00000       0.00000         0.00000       0.42410
## 378       0.00000       0.00000         0.00000       0.19010
## 379       0.00000       0.00000         0.00000       0.20510
## 380       0.00000       0.00000         0.00000       0.62410
## 381       0.03240       0.00360         0.03600       0.71650
## 382       0.00000       0.00000         0.00000       0.40280
## 383       0.00000       0.00000         0.00000       0.00000
## 384       0.00000       0.00000         0.00000       0.27550
## 385       0.00000       0.00000         0.00000       0.00000
## 386       0.00000       0.00000         0.00000       0.10940
## 387       0.00000       0.00000         0.00000       0.21560
## 388       0.00000       0.00000         0.00000       0.12820
## 389       0.01680       0.00450         0.02130       0.00000
## 390       0.00000       0.00000         0.00000       0.00000
## 391       0.00000       0.00000         0.00000       0.11340
## 392       0.00000       0.00000         0.00000       0.00000
## 393       0.03810       0.00260         0.04070       0.06780
## 394       0.00000       0.00000         0.00000       0.30480
## 395       0.02720       0.00320         0.03040       0.00000
## 396       0.00000       0.00000         0.00000       0.14127
## 397       0.00000       0.00000         0.00000       0.01650
## 398       0.00000       0.00000         0.00000       0.02260
## 399       0.00000       0.00000         0.00000       0.16820
## 400       0.03010       0.00090         0.03100       0.03500
## 401       0.00000       0.00000         0.00000       0.15950
## 402       0.00000       0.00000         0.00000       0.54120
## 403       0.00000       0.00000         0.00000       0.00000
## 404       0.00000       0.00000         0.00000       0.57490
## 405       0.00000       0.00000         0.00000       0.00000
## 406       0.00000       0.00000         0.00000       0.15560
## 407       0.00000       0.00000         0.00000       0.25330
## 408       0.01210       0.00140         0.01350       0.32090
## 409       0.00000       0.00000         0.00000       0.00000
## 410       0.00000       0.00000         0.00000       0.00000
## 411       0.00000       0.00000         0.00000       0.07500
## 412       0.00000       0.00000         0.00000       0.39750
## 413       0.00000       0.00000         0.00000       0.13500
## 414       0.00000       0.00000         0.00000       0.97110
## 415       0.00000       0.00000         0.00000       0.00000
## 416       0.00000       0.00000         0.00000       0.44120
## 417       0.00000       0.00000         0.00000       0.38860
## 418       0.00000       0.00000         0.00000       0.08530
## 419       0.06190       0.00430         0.06620       0.09690
## 420       0.00000       0.00000         0.00000       0.03760
## 421       0.03220       0.00900         0.04120       0.16740
## 422       0.00000       0.00000         0.00000       0.05730
## 423       0.00000       0.00000         0.00000       0.12560
## 424       0.00000       0.00000         0.00000       0.23380
## 425       0.00000       0.00000         0.00000       0.02970
## 426       0.00000       0.00000         0.00000       0.62390
## 427       0.00000       0.00000         0.00000       0.06630
## 428       0.00000       0.00000         0.00000       0.00000
## 429       0.00000       0.00000         0.00000       0.09820
## 430       0.00000       0.00000         0.00000       0.26320
## 431       0.00000       0.00000         0.00000       0.02860
## 432       0.00000       0.00000         0.00000       0.13500
## 433       0.00000       0.00000         0.00000       0.07830
## 434       0.00000       0.00000         0.00000       0.16210
## 435       0.00000       0.00000         0.00000       0.00000
## 436       0.00000       0.00000         0.00000       0.00000
## 437       0.00000       0.00000         0.00000       0.00000
## 438       0.00000       0.00000         0.00000       0.00000
## 439       0.00000       0.00000         0.00000       0.23260
## 440       0.00000       0.00000         0.00000       0.09420
## 441       0.02680       0.01160         0.03840       0.24500
## 442       0.00000       0.00000         0.00000       0.19780
## 443       0.00000       0.00000         0.00000       0.15270
## 444       0.00000       0.00000         0.00000       0.45260
## 445       0.00000       0.00000         0.00000       0.00000
## 446       0.00000       0.00000         0.00000       0.35680
## 447       0.00000       0.00000         0.00000       0.00000
## 448       0.02670       0.00160         0.02830       0.08550
## 449       0.00000       0.00000         0.00000       0.15620
## 450       0.00000       0.00000         0.00000       0.27920
## 451       0.02160       0.00660         0.02820       0.16840
## 452       0.00000       0.00000         0.00000       0.49780
## 453       0.00000       0.00000         0.00000       0.22620
## 454       0.00000       0.00000         0.00000       0.20090
## 455       0.00000       0.00000         0.00000       0.00000
## 456       0.00000       0.00000         0.00000       0.00000
## 457       0.02930       0.00660         0.03590       0.14020
## 458       0.02620       0.00370         0.02990       0.00000
## 459       0.00000       0.00000         0.00000       0.47020
## 460       0.00000       0.00000         0.00000       0.23670
## 461       0.01340       0.00180         0.01520       0.60500
## 462       0.00000       0.00000         0.00000       0.00000
## 463       0.04330       0.00300         0.04630       0.00000
## 464       0.00000       0.00000         0.00000       0.14650
## 465       0.00000       0.00000         0.00000       0.32680
## 466       0.02660       0.00420         0.03080       0.10670
## 467       0.00910       0.00190         0.01100       0.07940
## 468       0.00000       0.00000         0.00000       0.13500
## 469       0.00000       0.00000         0.00000       0.10830
## 470       0.00000       0.00000         0.00000       0.12260
## 471       0.01810       0.00130         0.01940       0.17840
## 472       0.00000       0.00000         0.00000       0.10020
## 473       0.00000       0.00000         0.00000       0.40660
## 474       0.00000       0.00000         0.00000       0.09090
## 475       0.00000       0.00000         0.00000       0.00000
## 476       0.00000       0.00000         0.00000       0.00000
## 477       0.00000       0.00000         0.00000       0.14550
## 478       0.00000       0.00000         0.00000       0.00000
## 479       0.00000       0.00000         0.00000       0.06810
## 480       0.00000       0.00000         0.00000       0.00000
## 481       0.00000       0.00000         0.00000       0.31500
## 482       0.00000       0.00000         0.00000       0.00000
## 483       0.00000       0.00000         0.00000       0.00000
## 484       0.00000       0.00000         0.00000       0.27620
## 485       0.00000       0.00000         0.00000       0.86750
## 486       0.00000       0.00000         0.00000       0.15160
## 487       0.00000       0.00000         0.00000       0.14930
## 488       0.00000       0.00000         0.00000       0.05990
## 489       0.00000       0.00000         0.00000       0.24720
## 490       0.00580       0.00290         0.00870       0.66190
## 491       0.00000       0.00000         0.00000       0.08700
## 492       0.00000       0.00000         0.00000       0.00000
## 493       0.00000       0.00000         0.00000       0.37210
## 494       0.00000       0.00000         0.00000       0.46990
## 495       0.00000       0.00000         0.00000       0.00000
## 496       0.03770       0.00560         0.04330       0.00000
## 497       0.00000       0.00000         0.00000       0.55170
## 498       0.00000       0.00000         0.00000       0.00000
## 499       0.06350       0.01650         0.08000       0.00000
## 500       0.00000       0.00000         0.00000       0.42340
## 501       0.03100       0.00510         0.03610       0.40080
## 502       0.00000       0.00000         0.00000       0.75670
## 503       0.00000       0.00000         0.00000       0.67150
## 504       0.03980       0.00500         0.04480       0.00000
## 505       0.00000       0.00000         0.00000       0.00000
## 506       0.01140       0.00120         0.01260       0.17340
## 507       0.02990       0.00720         0.03710       1.47650
## 508       0.04280       0.00660         0.04940       0.00000
## 509       0.00000       0.00000         0.00000       0.65230
## 510       0.03820       0.00470         0.04290       2.21290
## 511       0.00000       0.00000         0.00000       0.00000
## 512       0.00000       0.00000         0.00000       1.98050
## 513       0.00800       0.00000         0.00800       1.10640
## 514       0.05280       0.01170         0.06450       0.00000
## 515       0.00000       0.00000         0.00000       1.47190
## 516       0.05730       0.01580         0.07310       0.00000
## 517       0.05860       0.01120         0.06980       1.27780
## 518       0.00000       0.00000         0.00000       0.00000
## 519       0.00000       0.00000         0.00000       0.79800
## 520       0.00000       0.00000         0.00000       1.53660
## 521       0.00000       0.00000         0.00000       1.58920
## 522       0.01460       0.00180         0.01640       0.12930
## 523       0.00000       0.00000         0.00000       0.59390
## 524       0.02780       0.00640         0.03420       0.00000
## 525       0.00000       0.00000         0.00000       0.35320
## 526       0.00000       0.00000         0.00000       0.43340
## 527       0.00000       0.00000         0.00000       0.83420
## 528       0.00000       0.00000         0.00000       0.56180
## 529       0.00000       0.00000         0.00000       0.95550
## 530       0.01380       0.00290         0.01670       1.82680
## 531       0.00000       0.00000         0.00000       0.00000
## 532       0.07780       0.00800         0.08580       0.00000
## 533       0.00000       0.00000         0.00000       1.78700
## 534       0.02060       0.00460         0.02520       0.00000
## 535       0.00000       0.00000         0.00000       1.77900
## 536       0.00000       0.00000         0.00000       0.00000
## 537       0.00000       0.00000         0.00000       0.62810
## 538       0.00000       0.00000         0.00000       1.32370
## 539       0.05920       0.01120         0.07040       0.67900
## 540       0.00000       0.00000         0.00000       1.20680
## 541       0.00000       0.00000         0.00000       0.59440
## 542       0.00000       0.00000         0.00000       0.62600
## 543       0.00000       0.00000         0.00000       0.00000
## 544       0.02640       0.00470         0.03110       1.97150
## 545       0.01040       0.00170         0.01210       2.69650
## 546       0.02750       0.00410         0.03160       3.51890
## 547       0.00000       0.00000         0.00000       0.76370
## 548       0.00000       0.00000         0.00000       1.24970
## 549       0.04420       0.01040         0.05460       1.56030
## 550       0.00000       0.00000         0.00000       0.71700
## 551       0.00190       0.00160         0.00350       0.43900
## 552       0.00000       0.00000         0.00000       1.30790
## 553       0.05450       0.00520         0.05970       0.43160
## 554       0.00000       0.00000         0.00000       0.57820
## 555       0.03570       0.01080         0.04650       0.56500
## 556       0.00000       0.00000         0.00000       0.86740
## 557       0.00000       0.00000         0.00000       1.21300
## 558       0.00000       0.00000         0.00000       0.31220
## 559       0.07600       0.01040         0.08640       0.36080
## 560       0.00000       0.00000         0.00000       0.23680
## 561       0.03090       0.00820         0.03910       0.66100
## 562       0.00000       0.00000         0.00000       0.79290
## 563       0.00000       0.00000         0.00000       0.00000
## 564       0.00000       0.00000         0.00000       0.43580
## 565       0.00000       0.00000         0.00000       0.50580
## 566       0.00000       0.00000         0.00000       1.88470
## 567       0.03820       0.00720         0.04540       0.09400
## 568       0.00000       0.00000         0.00000       1.58130
## 569       0.00000       0.00000         0.00000       0.31140
## 570       0.00000       0.00000         0.00000       0.63910
## 571       0.00000       0.00000         0.00000       0.34870
## 572       0.00000       0.00000         0.00000       2.07340
## 573       0.04080       0.00470         0.04550       0.10910
## 574       0.04210       0.00700         0.04910       0.19380
## 575       0.00000       0.00000         0.00000       0.45160
## 576       0.04720       0.00840         0.05560       0.28540
## 577       0.00000       0.00000         0.00000       1.08290
## 578       0.01440       0.00350         0.01790       0.73870
## 579       0.00000       0.00000         0.00000       0.28350
## 580       0.00000       0.00000         0.00000       0.70130
## 581       0.00000       0.00000         0.00000       1.03630
## 582       0.04300       0.00930         0.05230       0.37020
## 583       0.00000       0.00000         0.00000       0.00000
## 584       0.00000       0.00000         0.00000       0.65670
## 585       0.00880       0.00090         0.00970       0.00000
## 586       0.02460       0.00330         0.02790       0.75260
## 587       0.01580       0.00540         0.02120       1.71170
## 588       0.00000       0.00000         0.00000       1.48570
## 589       0.04780       0.00630         0.05410       0.00000
## 590       0.00000       0.00000         0.00000       0.00000
## 591       0.00000       0.00000         0.00000       0.36170
## 592       0.00000       0.00000         0.00000       2.13290
## 593       0.00000       0.00000         0.00000       0.66780
## 594       0.02510       0.00460         0.02970       3.02780
## 595       0.01290       0.00640         0.01930       0.00000
## 596       0.00000       0.00000         0.00000       1.79770
## 597       0.00000       0.00000         0.00000       2.72430
## 598       0.00000       0.00000         0.00000       1.47240
## 599       0.00000       0.00000         0.00000       0.13030
## 600       0.00000       0.00000         0.00000       0.05780
## 601       0.03950       0.00840         0.04790       0.71150
## 602       0.01990       0.00290         0.02280       0.71100
## 603       0.00000       0.00000         0.00000       0.00000
## 604       0.00000       0.00000         0.00000       0.96310
## 605       0.04240       0.00850         0.05090       0.00000
## 606       0.00000       0.00000         0.00000       0.11960
## 607       0.00000       0.00000         0.00000       1.13270
## 608       0.00000       0.00000         0.00000       0.62850
## 609       0.04270       0.01200         0.05470       0.00000
## 610       0.00000       0.00000         0.00000       0.81020
## 611       0.00000       0.00000         0.00000       0.14460
## 612       0.00000       0.00000         0.00000       0.00000
## 613       0.04710       0.01050         0.05760       0.04700
## 614       0.00000       0.00000         0.00000       2.47650
## 615       0.00000       0.00000         0.00000       0.00000
## 616       0.00000       0.00000         0.00000       3.03630
## 617       0.03960       0.00000         0.03960       0.69670
## 618       0.00000       0.00000         0.00000       0.29100
## 619       0.00000       0.00000         0.00000       0.00000
## 620       0.00000       0.00000         0.00000       0.28710
## 621       0.00000       0.00000         0.00000       0.16570
## 622       0.01620       0.00170         0.01790       0.18540
## 623       0.00000       0.00000         0.00000       0.50660
## 624       0.00000       0.00000         0.00000       1.86340
## 625       0.00000       0.00000         0.00000       0.00000
## 626       0.00000       0.00000         0.00000       3.29120
## 627       0.04830       0.00680         0.05510       1.03610
## 628       0.03290       0.01400         0.04690       0.00000
## 629       0.00000       0.00000         0.00000       1.07190
## 630       0.07340       0.01560         0.08900       3.15560
## 631       0.00000       0.00000         0.00000       0.00000
## 632       0.00000       0.00000         0.00000       1.08420
## 633       0.00000       0.00000         0.00000       1.55460
## 634       0.00000       0.00000         0.00000       1.55010
## 635       0.04960       0.01200         0.06160       0.00000
## 636       0.00000       0.00000         0.00000       0.00000
## 637       0.00000       0.00000         0.00000       0.77010
## 638       0.00000       0.00000         0.00000       0.00000
## 639       0.00000       0.00000         0.00000       1.85760
## 640       0.00000       0.00000         0.00000       0.13240
## 641       0.02110       0.00410         0.02520       0.22330
## 642       0.00000       0.00000         0.00000       1.50170
## 643       0.02830       0.00530         0.03360       0.69590
## 644       0.06650       0.01300         0.07950       1.20800
## 645       0.00000       0.00000         0.00000       0.00000
## 646       0.00000       0.00000         0.00000       0.97640
## 647       0.00000       0.00000         0.00000       0.30640
## 648       0.00434       0.00300         0.00734       0.10810
## 649       0.00000       0.00000         0.00000       0.50970
## 650       0.03110       0.00510         0.03620       1.40070
## 651       0.02830       0.01030         0.03860       0.15870
## 652       0.03600       0.00247         0.03847       0.21570
## 653       0.00000       0.00000         0.00000       0.72530
## 654       0.00000       0.00000         0.00000       1.47360
## 655       0.00000       0.00000         0.00000       0.00000
## 656       0.09740       0.01020         0.10760       0.00000
## 657       0.00000       0.00000         0.00000       0.16810
## 658       0.05350       0.00840         0.06190       0.00000
## 659       0.00000       0.00000         0.00000       2.13070
## 660       0.07690       0.00780         0.08470       0.33380
## 661       0.00000       0.00000         0.00000       0.97320
## 662       0.00000       0.00000         0.00000       0.00000
## 663       0.00000       0.00000         0.00000       0.00000
## 664       0.00000       0.00000         0.00000       0.72560
## 665       0.00000       0.00000         0.00000       3.73300
## 666       0.01460       0.00110         0.01570       1.43600
## 667       0.01500       0.00510         0.02010       0.33210
## 668       0.00000       0.00000         0.00000       0.70570
## 669       0.00000       0.00000         0.00000       1.36380
## 670       0.00000       0.00000         0.00000       2.13170
## 671       0.03710       0.00000         0.03710       0.51460
## 672       0.00000       0.00000         0.00000       0.00000
## 673       0.07280       0.01340         0.08620       2.23580
## 674       0.00000       0.00000         0.00000       1.02030
## 675       0.05530       0.00940         0.06470       0.00000
## 676       0.00000       0.00000         0.00000       0.00000
## 677       0.08420       0.01230         0.09650       2.82360
## 678       0.02850       0.00460         0.03310       0.00000
## 679       0.03330       0.00760         0.04090       0.43680
## 680       0.00000       0.00000         0.00000       0.13980
## 681       0.00000       0.00000         0.00000       1.27660
## 682       0.00000       0.00000         0.00000       0.00000
## 683       0.00000       0.00000         0.00000       0.00000
## 684       0.00000       0.00000         0.00000       0.15530
## 685       0.00000       0.00000         0.00000       1.52910
## 686       0.00000       0.00000         0.00000       0.40020
## 687       0.03830       0.00000         0.03830       0.08470
## 688       0.00000       0.00000         0.00000       0.23920
## 689       0.00000       0.00000         0.00000       1.07630
## 690       0.00000       0.00000         0.00000       1.25280
## 691       0.00000       0.00000         0.00000       0.12860
## 692       0.00000       0.00000         0.00000       0.00000
## 693       0.05940       0.00390         0.06330       0.48940
## 694       0.00000       0.00000         0.00000       1.04530
## 695       0.00000       0.00000         0.00000       0.00000
## 696       0.00000       0.00000         0.00000       0.08600
## 697       0.05310       0.00380         0.05690       2.10530
## 698       0.05120       0.00950         0.06070       0.00000
## 699       0.00000       0.00000         0.00000       0.65810
## 700       0.02370       0.00400         0.02770       0.27220
recruit[is.na(recruit)] <- 0

recruit[,"eph.rate"] <- recruit$ephedra.emergence/3*100

## Percent of ephedra seedlings becoming established
mean(recruit$eph.rate)
## [1] 9.857143
## Percent of brome seedlings becoming established
mean(recruit$brome.prop.begin)
## [1] 51.36429